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SUMMARY 

With the coBopletion of an extensive study of global ozone patterns, 
this project has moved into the diagnostic phase of investigation of 
interhemispheric circulation differences. Early results are most 
interesting and will soon be joined by studies of the atmospheric mass 
structure inferred from the satellite measurements. 

The report includes a complete analyses of 0^ data inferred from 
NIMBUS-3 measurements, a discussion of future areas of study, description 
of the regression and inversion methods developed to infer atmospheric 
temperature and tropopause characteristics, as well as the plan to 
process the satellite data for a systematic study of the relative 
circulation differences between northern and southern hemispheres. 


i 
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1.0 INTRODUCTION TO GRANT OBJECTIVES 

It is the purpose of this research project to; 

a) Investigate, with the aid of satellite measurements, the 
fundamental processes and interactions affecting inter- 
hemispheric differences in mass structure and energetics, and 

b) To apply measurements from satellites to the problem of 
large-scale global weather forecasting. 

In this, the second, annual report we present results of studies 
using measurements of infrared radiance from the Infrared Interferometer 
Spectrometer (IRIS) flown on NASA's NIMBUS-3 satellite. These data, 
with the aid of NASA, have been analysed selectively to study the x , 
y , t distribution of atmospheric ozone and its relation to large- 
scale circulation parameters. The satellite data provide a global 
view during the period April-July, 1969. Dr. James Lovill, as his 
dissertation, has examined this ozone climatology, compared it with 
ground-based observations, and used it to infer the characteristics of 
regional and planetary scale circulation systems in both the Northern 
and Southern Hemispheres. In addition to the study of mid-latitude 
waves, reported herein, 0^ data are used to infer the magnitude of 
jet stream wind maxima and for other special studies. 

From the same atmospheric spectra used to derive 0^ , selected 

spectral radiance measurements have been used by Mr. S. Srivatsangam to 
obtain vertical temperature profiles in both hemispheres from a statistical 
inversion method developed by Rodgers (1969) . Together with a less complex 
regression method to determine upper tropospheric temperature from 
radiance, the satellite results are "calibrated" over the Northern 
Hemisphere using standard NMC data. Examples of the results are included 
in this report. 
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Mr. Robert Adler has devised a special technique to detect tropo- 
pause temperature and pressure using IRIS data. This allows rapid infer- 
ence of inter-hemispheric circulation differences and may be a great aid 
to operational use of satellite data. Both he and Srivatsangam plan to 
use the extensive satellite data set to study the geopotential height 
fields, the available potential and kinetic energy magnitudes of each 
hemisphere. 

Special data used in support of the satellite measurements include 
ground-based ozone measurements, conventional meteorological analyses over 
the Northern Hemisphere and a 5-year set of upper air wind data from the 
Australian sector. Further work will include ^elective ^hopper Radiometer 
(SCR), Satellite InfraRed ^ectrometer (SIRS), and IRIS measurements 
from the NIMBUS-4 as well as results from the HOLE experiment. 

The following sections provide a more detailed discussion of grant 
research. It is without doubt that a controlled use of the recently 
available satellite measurements to study the circulation differences 
between Northern and Southern Hemispheres offers atmospheric scientists a 
potential quantum-jump in knowledge regarding the significant processes 
and interactions and their natural variation. 



4 


2.0 DISCUSSION OF RESULTS 

2.1 Use of Ozone to Study Atmospheric Transport Processes 

2.11 Meridional 0^ ^adient vs. jet stream maxima; temporal and 
geographical variations of ozone 

Studies in recent years have shown distinct differences of the 
atmospheric general circulation of the Northern and Southern Hemispheres. 
Earlier publications postulated a more zonal flow in the Southern 
Hemisphere supposedly because of less land mass and fewer mountain 
massifs. More recent investigations show that such statements need to 
be qualified. 

Wooldridge and Reiter (1970) have shown that significantly 
stronger horizontal anisotropy of flow prevails at cyclone wavelengths 
in the Southern Hemisphere than in the Northern Hemisphere. In the 
Southern Hemisphere meridional perturbations exceed zonal flow perturba- 
tions. Eulerian spectral densities of the zonal component in the 
Southern Hemisphere were seen to be only half of those in the Northern 
Hemisphere. This would seem to imply that there is less zonal energy 
at cyclone-wave frequencies in the Southern than in the Northern 
Hemisphere. This may occur because of a lack of orographically induced 
planetary longwaves in the Southern Hemisphere that are superimposed 
upon the hemispheric jet stream zones. Apparently the relatively 
zonal character of the Southern Hemisphere pertains to the long 
planetary-scale waves only. 

Interest in global total ozone measurements is widespread and 
increasing. The paucity of total ozone measuring stations over the 
oceans and in the tropics in the Northern Hemisphere and throughout 
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the Southern Hemisphere has greatly hindered a complete understanding 
of the role that stratospheric ozone plays in the general circulation 
processes of the atmosphere. With total ozone sensors (such as the 
IRIS instrument) on earth-orbiting satellites the ozone over the 
entire earth is measured daily. The conclusions of some of the 
analysis from the first satellite (Nimbus III) capable of these 
measurements follow. 

A statistical analysis of the total ozone data for two days 
indicated that average standard error between data obtained by the 
inversion technique and data obtained by regression coefficients was 
only 3.7 percent. When surface total ozone measurements were compared 
with regression data at monthly intervals a standard error of 3.6 
percent was seen. 

A strong correlation was discovered between the meridional gradient 
of total ozone and the wind velocity in jet stream systems. A graph 
has been constructed which will allow the computation of the wind 
velocity at 250 mb if the total ozone distribution in the vicinity is 
known. This should be a very useful relation to apply over the oceanic 
regions of the Northern Hemisphere and throughout the entire Southern 
Hemisphere. Wind velocities at jet aircraft operation levels could 
then be forecasted by consulting the daily satellite total ozone analysis. 
This method to obtain wind velocity might ideally complement the technique 
to infer upper troposphere wind direction, based on water vapor measure- 
ments, being developed at NASA. 

A study of the total ozone distribution over a Pacific Ocean extra- 
tropical cyclone and its associated jet stream indicated large vertical 
motions and strong horizontal advective processes were present near this 
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system. Large amounts of total ozone were measured immediately behind 
the system and low amounts in front which were indicative of air rich in 
ozone moving southward and downward behind the system and air poor in 
ozone moving northward and upward in front of the storm. 

Time-longitude variations of the total ozone for several latitudes 
in the Northern and Southern Hemisphere for June 1969 were analyzed. 
Several interesting facts were ween in these results. Firstly, one can 
observe that the speed at which the 'waves' progress eastward is greater 
in the winter hemisphere (here the Southern Hemisphere) . Secondly, the 
speed of eastward progression decreases as one approaches the lower 
latitudes in the winter hemisphere. Thirdly, in the equatorial region 
and in the Northern Hemisphere there is not a strong eastward progres- 
sion of 'waves' but at best only weak eastward movement. Progression 
rates as high as fifteen degrees of longitude per day were seen at the 
high latitudes in the Southern Hemisphere. As one follows the 'waves' 
in time, continuous changes in intensity are observed (see also 2.12). 

When total ozone data for the period of observations are plotted 
against latitude for the various months the following is indicated. 

For the first time the average distribution of ozone has been described 
during the automn and winter at close intervals, both temporally and 
spatially, in the Southern Hemisphere. The smallest variation of 
total ozone from spring to summer was seen at 15°N. The lowest mean 
total ozone value for the globe was at 6°S. 

When two maps were prepared describing the distribution of total 
ozone for the entire period for the Northern and Southern Hemispheres, 
the following was seen. Closely spaced total ozone isolines appear to 
have revealed a reasonably strong jet stream west of South America at 
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approximately 35°S. This system, at least at the intensity indicated 
by the total ozone gradient, has gone more-or-less unnoticed in the 
past. 

Additionally, a large region of low total ozone was observed in 
the tropics, extending from southeast Asia eastward to the central 
Pacific. Earlier studies have indicated that this region also has an 
anomalously high albedo. The manner by which this correlates with low 
total ozone concentrations is not at all clear. 

Both the Northern and Southern Hemispheres appear to provide 'an- 
choring' mechanisms for total ozone ridges. The most predominant of 
the ridges is also usually associated with a 'tight' ozone gradient. 

In the Northern Hemisphere these locations are eastern North America, 
central Europe, western Asia, and eastern Asia. The ridge over the 
Soviet Union appears to be the weakest of the four ridges. In the 
Southern Hemisphere 'ozone' ridges are evident over the east Indian 
Ocean, the central Pacific Ocean and the eastern Atlantic. In general, 
the ridge pattern is easier to delineate in the Southern Hemisphere than 
in the Northern Hemisphere. 

The 'tightest' gradient of total ozone in either hemisphere was 
seen over Japan. This correlates closely with the high frequency of 
jet streams observed over the region. 

An average global value of the total ozone for the entire period 
was calculated to be 330.3 m*atm*cm. 

A paper containing these results in detail is enclosed as Appendix 
I to this report. 
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2.12 Extended study of mid-latitude wave patterns 

Special time-space displays of the ozone data are being used in 
an extended study of mid-latitude wave patterns. This work will be 
summarized in a special report, now in preparation. 

2.2 Atmospheric Thermal Structure From IRIS Spectra 

Recent research under this grant (Reiter, et al., 1971) using 
some 88 near-simultaneous (At^ 3 hours) radiosonde and IRIS radiance 
observations indicate that in the stratosphere the correlations between 
the ambient temperatures and the radiances of the infrared wave-numbers 
weighted predominantly in this atmospheric layer (such as 675 cm ^) are 
so significant that temperatures inferred directly from the scatter- 
diagrams are comparable in accuracy to radiosonde observations. Differences 
between the inferred temperatures and the NMC grid temperatures a few 
hours apart are of the order of 2°C. 

However, isobaric level temperatures in the troposphere and the 
radiances of the IR wave-mombers weighted significantly about these 
levels show such a wide scatter that direct inference of temperatures 
from the scatter-diagrams is impossible. This difficulty leads one 
to the use of other, and generally computationally more costly, methods 
of analysis in order to study stratosphere-troposphere interaction using 
satellite data. 

The method due to Rodgers (1970), which we have decided to employ, 
uses estimation theory as a means to invert the linearised radiative 
transfer equation to obtain atmospheric vertical temperature profiles. 

Two important constraints are imposed by this method on the resulting 
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data for analysis. One is its limitation to clear-sky cases only and 
the other is the satisfaction of the linearity criterion mentioned 
above. The differentiation between cloudy and clear sky cases is 
discussed in both 2.21 and 2.23. 

More important than the vertical thermal structure in the analysis 
of tropospheric circulation features is the distribution of mass. 
Despite the inter-relationship between the mass and temperature fields 
"there is no method of inferring the sea- level pressure field hydro- 
statically from the three-dimensional temperature field alone" (Lorenz, 
1967, p. 28) and, because of this restriction it would not be possible 
to obtain the geopotential height of a given isobaric level from the 
temperature profile only. Statistical analysis (Smith and Fritz, 1969) 
shows that this is not an unavoidable condition. 

In our own analysis here at Colorado State University we intend 
to derive, at the present time, only the relative height fields 
(Z' -fields) of some 5 isobaric levels from 100 mb to 700 mb. The Z'- 
fields (deviations from the mean Z for an area S comparable to the 
area of a hemisphere) are expressible in terms of fast-converging 
empirical orthogonal functions Fj^ and their coefficients Zj^ as 

n 

Z' (x, y, p, t) = I Z (x, y, t) F (p, t) (1) 

k - 1 

The analytical solution of equation (1) as given by Holmstrom (1963; 
see also Obukhov, 1960) is 
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We have developed a program for evaluating and Fj^ , and used it 

for computations in the range k ^ 6, P = 1000 mb to 100 mb and P = 700 
mb to 100 mb. 

By using the hydrostatic equation and the equation of state in (1) 
we obtain 


f = -z-iigi-Fk 

^ ^ R dp 

(3a) 

= ^ Z, E, 

, k k 

(3b) 


where are defined through (3a) and t' are defined as temperature 

deviations from the area mean, analogous to Z' (see Smith et al., 1970 
and Holmstrom, 1963). From (3b) and (1) the Z'-field is obtained from 
the t' -field through the matrix equation 

Z* = t' E'^ F (4) 

We have solved equation 4 for the data of 10 July 1969. The first 
results indicate considerable deviations from actual values. It is 
not clear whether the 5-mode solution is over-descriptive or under- 
descriptive. Holmstrom (1963) says the functions Tj^ (p) which he 
develops for the temperature-field from Fj^ (p) are not fast-converging. 
On the other hand Z^'s are very small, especially when Fj^ are nor- 
malized in the interval 700 mb-100 mb. We are trying to evaluate 
equation (4) with a 4-mode system. 

If subsequent results are encourageing we intend to use empirical 
orthogonal functions to the vorticity tendency equation, first with 
conventional data and later with IRIS-derived Z'-field data. Initial 
theoretical considerations indicate that such an analysis would enable 
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the separation of different atmospheric developmental processes such 
as the development of the semi -permanent high and lows and baro- 
clinic development due to ultra-long and cyclone waves. 

Should the empirical orthogonal function approach prove not to be 
a feasible one, our circulation studies would still be based on: 

a) studies of the time and space variation of the mean temperature 
of atmospheric layers , and/or 

b) geopotential height fields derived using a regression method 
very similar to the one employed operationally with the SIRS 
data 

Throughout the remainder of this section we report on results of 
inferring the atmospheric temperature at distinct levels . This is a 
fundamental step that allows us to compare with other work of this type. 
However, for our circulation study, the comparisons with NMC data must 
be regarded as worst-case since random errors in temperature retrieval 
that affect RMS level comparison are removed from atmospheric layer 
(mass) comparisons. 

2.21 Results from radiometric inversion and from direct regression on 
radiance values 

The InfraRed Interferometer spectrometer, IRIS, is a Michelson 
interferometer employing a beamsplitter which divides the incoming 
radiation into two approximately equal components, one of which is 
made to reflect from a stationary mirror and the other from a moving 
mirror. The interference of the two beams produces an interference 
pattern (interferogram) from which a radiation spectrum is constructed 
by applying the inverse Fourier transform. Spectral radiances in 
wave-number bands of 5 cm ^ can be obtained from this spectrum. The 
weighting functions for these bands in a vertical atmospheric 
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coliunn are given by 


dx . 

K = 1 

i d iln p 


( 5 ) 


where is the transmissivity o£ the atmosphere in the i^^ wave- 
ntunber interval. The weighting functions for such wavenumber intervals 
in the 15ym CO2 band, some of which we used to retrieve atmospheric 
temperatures by Rodgers* (1970) method, are given in Figure 1. 

Our aim was to obtain temperatures at 6 isobaric levels, viz., 

700 mb, 500 mb, 400 mb, 200 mb, 100 mb and 30 mb. Cloud interference 
makes the retrieval of tropospheric temperatures difficult. One solu- 
tion to the problem is to prescribe sufficiently accurately, a set of 
parameters such as the surface temperature in the field of view of 
IRIS and cloud-top heights. (For other usable sets of parameters see 
Chahine (1970)). In an iterative procedure such as the one described 
by Smith et al. (1970) — where cloud-top heights are treated as variables — 
the only parameter needed, in addition to the spectral radiances them- 
selves, is the surface temperature. This can be accomplished over oceans 
by using independent satellite observations whereas, over land the 
diurnal temperature variation makes this rather difficult to carry out. 

Another way of solving the cloud problem would be to prescribe 
adequate lower limits for spectral radiances whose weighting curves for 
the gaseous atmosphere peak in the lower troposphere and therefore would 
be affected significantly by clouds at any level in the atmosphere. 

Thus clear-sky and cloudy areas in a region (such as a latitude band) 
could be distinguished. 745 cm"^ is a wave-number that satisfies this 
criterion (see Figure 1) . Radiances in the window region, which are 
influenced primarily by surface and cloud-top temperatures could also 
be used this way, but these depend very much on surface characteristics. 
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Fig. 1 - Atmospheric weighting functions for 667 cm (15ym) CO2 absorption 
band. The curves were calculated assuming an instrument function 
5 cm" wide at the half maximum point, centered on the wave number 
indicated. (After Conrath, et. al., 1970). 
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TABLE 1 


Limiting values of 745 cm ^ radiances for distinguishing between clear 
and cloudy areas. is equivalent black body temperature in degrees 
centigrade. 


Initial Values 


Values when lower tropospheric 
deviations are random 


Latitude 

Band 

(Degrees) 

745 cm ^ 

-1 -1 -1 
erg cm s sr 

T 

bb 

Degrees C 


745 cm ^ 

-1 -1 -1 
erg cm s sr 

T 

bb 

Degrees 

0-30 

80 

-14.0 

90 

- 6.6 

30-60 

78 

-15.6 

85 

-10.2 

60-81 

68 

-23.6 

80 

-14.0 
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Since the satellite observations are made (in the case of NIMBUS III 
which carried IRIS) very nearly at local noon and midnight except at 
subpolar latitudes allowance has to be made for the diurnal temperature 
variation over land. Wave-number 745 cm ^ is much less dependent upon 
surface thermal characteristics and obviates the need for ocean surface 
temperature and therefore was used. The limiting values used for 
different latitude bands are listed in Table 1. Atmospheric temperatures 
obtained under such criteria and using Rodgers' (1970) method were 
interpolated linearly from the subsatellite track to a 5° latitude by 
5° longitude grid for the Northern Hemisphere. The data were again 
interpolated using a 16 point Bessel interpolation to fill the NMC grid 
points. Contours fitted for these data using a CDC user program were 
plotted on microfilm by a CDC 6400 Con^uter. From these analyses, 
data for locations corresponding to 20 Eurasian and 18 North American 
radiosonde stations were collected and plotted against the actual 
radiosonde data for each of the levels mentioned above. (The North 
American radiosonde data are for 1200 GMT and the Eurasian ones for 
0000 GMT; so the local time is approximately the same for Asian and 
American data.) The scatters were studied for any consistent (bias) 
deviations of the IRIS-derived temperatures from the radiosonde 
observations. Since we started with low limits for 745 cm~^ radiances 
the initial IRIS-derived temperatures were colder than radiosonde 
temperatures at the 700 mb and 500 mb levels. So the limiting 745 cm"^ 
radiance values were raised, and the new scatters checked and the proce- 
dure repeated until the errors at the lower tropospheric levels were 
con 5 >letely random. Under such conditions the RMS deviation of the IRIS- 
derived temperatures is sufficiently independent of any cloud influence. 
The values of 745 cm ^ radiances for which the errors at the 700 mb 
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TABLE 2 


RMS and mean deviations for the data presented in Fig. 2, which are for 
days 3 and 8 July, 1969. 


Level 

(mb) 

RMS Deviations 
(Deg. C) 

Mean deviation of IRIS-derived 
temp, from radiosonde temp. 
(Deg. C)* 


Eurasian data 

All data 

Eurasian data 

All data 

700 

3.61 

3.78 

0.29 

0.11 

500 

2.76 

2.87 

0.18 

0.31 

400 

3.14 

3.46 

0.48 

0.99 

200 

4.18 

3.78 

-0.38 

-0.92 

100 

2.39 

2.94 

0.08 

-0.25 

30 

2.15 

3.24 

0.92 

2.16 


* 

Hence a positive value indicates a wanner radiosonde temperature compared 
to IRIS-derived temperature. 
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and 500 mb levels were con 5 )letely random are presented in Table 2. 

The r.m. s. error is much larger at the 700 mb level coD 5 )ared to the 
r.m.s. error at the 500 mb level. The final IRIS-derived ten^jeratures 
plotted against radiosonde ten^eratures for all the levels considered 
are presented as Figures 2a - 2f. 

The r.m.s. errors are generally smaller for the Eurasian data 
than for the combined Eurasian-North American data, as seen from Table 
2. Since there is no reason to expect that the analytical method 
employed here has any regional bias, this increase in error for the 
combined data might be attributed to the time difference between the 
satellite and radiosonde observations. This difference is larger for 
North America and smaller for Eurasia, but is close to zero even 
here only at 0000 GMT. Thus, at least part of the r.m.s. errors of 
Table 2 might be due to time difference in observations. 

In our earlier report (Reiter, et al., 1970) a scatter diagram 
for 100 mb temperatures against v = 675 cm ^ radiances was given. For 
that scatter a third degree curve seemed to be in best fit. 100 mb 
temperatures for 2 days in July 1969 were obtained from this curve. 
Distributions of these temperatures in the NMC grid area are presented 
here (Fig. 3a and 3e) and compared to the NMC temperature distribution 
(Fig. 3c, 3d, 3g and 3h ) at 0000 GMT and 1200 GMT on the same days. 
The distributions obtained by using Rodgers' (1970) method are also 
presented (Fig. 3b and 3f) . A con 5 >arison of these three types of 
distributions indicates that in general there is qualitative agreement 
between them both in temperature gradients (and hence in thermal wind 
shears) as well as in contour orientation (i.e., trough-ridge location). 
The results obtained are good considering the inability of the present 
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FIGURE 2a 
Level; 700mb 

Setter 0 ^ ]R!S-derived temp. (Deg. C) agoif'st Radiosonde temp. (Ocg. C) 
OS E'JRASIA:; PADIOSOTCES (0000 GMT) 

*+ ri.AMERlCAn RADlOSOrDES 0200 GMT) 




FIGURE 2b 


Level; 500mb 

Setter of IPiS-derived temp. (Deg. C) against Rodiosonde temp. (Deg. CJ 
90 E'JRASlArt RAD10S0M)ES (0000 &MT) 

!l.A^tERICAri RAD;0S0M)ES (!200 GMT) 




Radiosonde Timperotorc 
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FIGURE 2c 
Level; 400mb 

Satter of IRiS-derived temp, (Deo. C) ogoinst Radiosonde temp. (Deg. C) 
00 E'JRASUN RAD]0S0tX)S;S (0000 GMT) 
ri,A^f:RlCA^i RAD10S0fX)ES (1200 GMT) 




Radiosjcmde Temperatures 
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FIGURE 2d 


Level: 200mt» 

Satter of IR'iS-derived temp. (Deg. C) against Radiosonde temp. lOeg. C) 
00 £!jPASIArj RADIOSONDES (0000 GMT) 

4+ ^'.A^^;RlCArJ RADI0S0^DES (1200 GMT) 

- 50 , 


OO. 


-* 0 , 




Radiosonde Temperatures 
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FIGURE 2e 
Leve ! I 1 OOmb 

Setter of IRlS-derived temp. (De 9 - C) agaiost Radiosonde temp. (Oeg. C) 
00 EURASIAN radiosondes (0000 GMT) 

++ N. AMERICAN RADIOSONDES (1200 GMT) 

- 40 . 


-SO. 
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FIGURE 2f 
Level: 30mb 

Satter of IR I S-de r i ved temp. (Deg. C) ogoinst Radiosonde temp. (Oeg. C) 
00 ELRASlAfJ RAOIOSON)ES (0000 GMT) 

++ rLA^t:R!CA^i RAOlOSOfCES (1200 GMT) 
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Fig. 3a - 100 mb temperatures (Degrees C.) obtained by correlation between 
V = 675 cm radiances and 100 mb temperatures. 3 July 1060. 
Number of data points 1597. 
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Fig. 3c - Northern Hemisphere Radiosonde terrperatures (Degrees C.) at 
100 mb. National Meteorological Center (NMC) mapping. 

8 July 1969. 0000 GMT. 





Fig. 3d - Northern Hemisphere Radiosonde temperatures (Degrees C . ) 

100 mb. National Meteorological Centre (NMC) mapping. 

3 July 1969. 1200 GMT 
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Fig. 3e - 100 mb temperatureSj^ (Degrees C.) obtained by correlation 
between v = 675 cm radiances and 100 mb temperatures. 
8 July 1969. Number of data points 1193. 
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Fig. 3i - Mean July 1969 tropopauses along North American coastal regions. 

A tertiary region of stability exists at about the 140 mb level 
at most latitudes considered here. 
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methods for inverting the radiative transfer equation to detect 
temperature inversion in the upper atmosphere. Smith et al., (1970) 
discuss the problem; sec. 2.3 includes a suggested method of improve- 
ment in the tropopause region. The tropical tropopause which extends 
into very high latitudes during the summer (Fig. 3i) causes some of 
the deviations of the derived temperatures from the NMC temperatures. 

NMC temperature distributions indicate very weak temperature 
gradients over the oceans. Satellite observations also indicate such 
gradients over oceans, but the differences between the oceans and the 
continents are not as large as in the case of the NMC distributions. 

These discrepancies are probably due to the very meager radiosonde 
data available in oceanic regions. Thus, IRIS-derived temperatures 
show greater zonal symmetry than the NMC data. 

2.22 Special comparison of results with 30 mb radiosonde data 

Figure 2f shows that at the 30 mb level the North American 
radiosonde temperatures are relatively warmer compared to the IRIS- 
derived temperatures whereas the Eurasian data are randomly distributed 
about the one-to-one correspondence line. (See also Table 2: The 

r.m.s. error for the 30 mb level temperatures is much smaller for the 
Eurasian data compared to the combined Eurasian-North American data.) 

This could be due to a natural diurnal temperature variation or a 
larger radiational error of the instrument used in North America 
compared to that of the Eurasian instnunent. The former explanation 
is difficult to maintain in view of the fact that the Asian stations, 
which have made the observations at about the same local time as the 
North American stations, do not indicate such large systematic departures 
from the IRIS-derived temperatures. 
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One way to cheek any such large radiational errors would be to 
compare the radiosonde observations made at 0000 GMT and 1200 GMT 
along the two coastal regions of the North American Continent where 
the Local Standard Time (LST) for the above observational hours are 
significantly different. These LST equivalences are given in Table 3, 
after making an allowance of 2 hours for the ascent of the balloon to 
the 30 mb level. Table 4 contains the solar elevation angles in July 
at different latitudes along these coasts an hour after the release 
of the balloon; the radiational error generally accumulates after 
this time in the presence of sun because the balloon is in the strato- 
sphere exposed to high insolation is not well ventilated. We will 
be concerned only with the latitude band 30-50N. From Tables 3 and 4 
it is seen that the observations numbered 1 and 4 in Table 3 are made 
under low solar elevations and those numbered 2 and 3 under high solar 
elevations. 

The mean temperatures for July 1969 at 30 mb at 4 east coast 
and 4 west coast stations for 0000 GMT and 1200 GMT were calculated 
and are presented in Fig. 4a. By averaging data for an entire month 
we hope to have eliminated all the diurnal temperature variation due 
to changes of synoptic situations. In Fig. 4a all observed values are 
consistent with the fitted profiles except Pittsburgh at 0000 GMT. 

It is readily inferred from the diagram that the daily mean temperatures 
(in July 1969 at 30 mb) are very nearly the same in the two coastal 
areas . 

Except for some small deviations the observations made imder low 
solar elevations coincide in Fig. 4a and the observations made under 
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North American coastal regions. 
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TABLE 3 


Approximate local standard time (LST) for the 30 mb radiosonde obser- 
vations in North American continental coastal areas between 30-50N. 
Ascent to this level assumed to take place in 2 hours. Arbitrary ob- 
servation numbers are given in parentheses. 
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TABLE 4 

July solar elevation (to the nearest full degree) for 0000 GMT and 
1200 GMT radiosonde observations an hour after the release of balloons 
at the Local Standard Time (LST) given. Arbitrary observation numbers 
written in parantheses. W.C. - West coast of North America. E.C. - 
East coast of North America. Data from Smithsonian Meteorological Tables 
(R.J. List, editor). 


Latitude 

Degrees 

North 

0000 

1200 

GMT 

LST 

E .C. 

2000 

(1) 

W.C. 

1700 

(2) 

E.C. 

0800 

(3) 

W.C. 

0500 

(4) 

50 

0 

27 

36 

8 


45 

0 

26 

36 

6 


40 

0 

25 

36 

3 


35 

0 

24 

36 

1 


30 . 

0 

23 

36 

0 



NOTE: Refraction of solar radiation would advance sunrise and delay sun- 

set. This factor has not been taken into account in the above calculations. 
Because of this factor more radiation would fall on the radiosonde at 2000 
LST and 0500 LST than is suggested by its solar elevation angles given 
above (see Mitra, 1951). 





LOCAL STANDARD TIME (hours) 

Fig. 4b - Apparent diurnal variation of temperature (Degrees C.) at 30 mb at different latitudes along 
North American coastal regions in July 1969 as deduced from Fig. 4a. 
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high solar elevations coincide also. The latter are some 2C warmer 
than the former south of 40N and gradually decrease with latitude 
north of that latitude. 

Treating the two coasts as a single meridional sector in view of 
the invariance of the daily mean tenq>erature mentioned above the 
apparent diurnal variation of temperature at 30 mb as obtained from 
the radiosonde data of Figure 4a are plotted as a time-series in Figure 
4b. It is seen from this diagram that the diurnal maxima and minima 
must occur at 1400 LST and 0100 LST respectively at all latitudes 
between 30N and SON (to a reasonable degree of approximation) . This 
periodicity is extremely odd, both the increase in temperature from 
0100 LST to sunrise and the decrease from 1400 LST onwards being quite 
inexplicable. The amplitude of the apparent diurnal variation is in 
excess of 6C in the latitude band 33-38N and decreases to some 2.8C 
at 46N. 

The following remarks by Craig (1951) concerning the heating 
rates in the ozone layer and below are of interest: "Karandikar's (1946) 

computations show that the heating produced by infrared absorption of 
solar energy can be safely neglected above 30 km in comparison with the 
heating produced by ultraviolet absorption of ozone. Below 30 km, 
on the other hand, the former process becomes dominant, water vapor 
playing the most important role. . . The absolute values of the rates 
of heating and cooling are relatively low, of the order of magnitude of 
O.IC per three hours below 30 km and of 1C per three hours at 35-50 km." 
The radiosonde -given temperature change at the 30 mb (- 25 km) level, 
as seen in Fig. 4b, exceeds even the latter value in the 33-38N latitude 
band. 
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Finger and Mclnturff (1968), using North American radiosonde 
data for the period Jan. 1964 to June 1966 at the 30 mb level and above, 
have evaluated the diurnal temperature range R(T) considering stations 
which make both 0000 GKT and 1200 GMT observations either in the sun 
or in darkness. They arrive at these conclusions: The computed diurnal 
temperature range R(T) shows a definite increase with altitude in 
summer, R(T) at 10 mb being double that at 30 mb and that at 5 mb being 
2.5 times the R(T) at 10 mb. In winter the change with height is not 
so well-defined. At 30 mb R(T) is less than 0.5C in summer and generally 
decreases with latitude (Fig. 4c). The true diurnal change most 
probably peaks near sunset and has a minimum near sunrise. These 
conclusions are in very good agreement with Craig's (1951) remarks above. 

Fritz (1970) states that the diurnal variation of stratospheric 
temperature (integrated in the vertical) is small and is of the order 
of 0.5 erg cm ^ sr ^ s The temperature equivalent for this quantity 
is approximately 0.5C for a black-body. The integration referred to 
here is over the altitude range in which v = 669 cm ^ has a non-zero 
weighting function value. Our own studies using v = 675 cm ^ which has 
maximum weighting in slightly lower stratospheric levels (Fig. 1) 
show that there is very little diurnal variation of radiances in this 
wave-number. As this is probably also true of individual levels above 
the 100 mb level the differences between the 0000 GMT and 1200 GMT 
observations seen in Fig. 4a are most probably radiational errors in the 
sonde. 

Figure 5 also lends support to this contention. In this diagram 
the 1200 GMT east and west coast temperatures at 30 mb are differentiated 
and plotted against IRIS-derived temperatures which are a good approxima- 
tion to the daily mean value since the observations are made at local 
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Fig. 4c - Summertiine variation of diurnal temperature range 

R(T) in degrees C. for 30, 10 and 5 mb as a function 
of latitude. (From Finger and Mclnturff, 1968). 
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IRIS TEMPERATURE (*0 


Fig. 5 - North Americcin West Coast and East Coast radiosonde temperatures 
at 1200 GWT compared to IRIS-derived temperatures. Data for 
3, 8 July 1969. Local Standard Time (LSI) given is adjusted by 
2 hours for balloon ascent to the 30 mb level. 
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noon and midnight. It is seen that the east coast soundings (at 0900 
LST) are generally warmer than the west coast soundings (at 0600 LST) 
even on individual days. 

If these rather large apparent diurnal variations are adjusted, 
the IRIS-derived and radiosonde temperatures at the 30 mb level would 
agree very well over all the regions where comparisons have been made. 

2.23 Additional comparison of temperature -radiance relation for cloud- 
contaminated and "cloudless" data 

Two sets of "simultaneous" IRIS radiance and NMC temperatures were 
compiled. One set of data consisted of data from which cloud- contaminated 
points had been eliminated by a technique to be described. The other 
set retained data from points affected by clouds. 

The method of removing cloud-affected data is as follows. For each 
day of data (there were three days of data in each set) frequency distribu- 
tions of radiance in the 735 cm ^ channel, N(735), are made for each 
ten degrees of latitude from 20°N to 80°N. The median of each of the 
frequency distributions is calculated. From comparisons of the spatial 
variations of N(735) and video cloud information, a criterion for elimina- 
tion of cloud-affected data was subjectively determined. The criterion 
eliminates data points where N(735) is less than N(735) - 8.0, where 
N(735) is the median for the latitude band of the data. This criterion 
was applied for each of three days of data (3, 4, 8 July 1969) to obtain 
the "cloudless" data set. The data were averaged over 5 points. The 
other set, the"with- cloud" set, was obtained by setting the cloud elimina- 
tion criterion low enough so that cloud-affected data is accepted into 


the data set. 
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The temperature- radiance correlation matrix for the "cloud-free" 
set and the cloud-affected data set are given in Table 5. The "cloud- 
free" correlations are the upper ones in each box; the ’tloud-affected" 
correlations are the lower ones. The differences between the"cloud-free" 

and "cloud-affected" correlations are less than 0.01 for channels 665 cm 

-1 -1 -1 
to 695 cm . For channels 705 cm to 900 cm the differences vary from 

0,03 to 0.45 . This relation is to be expected from the distribution 
of the individual chanel weighting functions. The channels with large 
differences have weighting functions with peaks in the middle and low 
troposphere. Those are the channels that are most affected by the pre- 
sence of clouds. The "cloud-affected" set had 567 data points, the 
"cloud-free" set, 469. Thus about twenty percent of the data were 
eliminated by judging it cloud-contaminated. 

The correlation matrix for the"cloud-free" data displays the ex- 
pected high positive correlation between low wave number radiance and 
high altitude temperature (e.g., N(675) and 100 mb temperature) and 
between high wave number radiance and low altitude temperature (e.g., 
N(725) and 850 mb temperature). The data set included measurements 
from 20°N to 80°N. Therefore, there is a strong negative correlation 
between tropospheric and stratospheric temperature. This is reason 
for the large negative correlations in the botton-left and upper-right 
comers of the correlation matrix. 

Regression equations for the two sets of data for temperatures at 
various pressure levels are given in Table 6. The fraction of the 
variance explained and the standard error are also given. A stepwise 
regression procedure was used with the selection of variables halted 
when the percent variance explained failed to increase by at least 1%. 



Temperature at Pressure Levels 


Temperature vs. Radiance Correlation Matrix 



665 

675 

685 

690 

695 

700 

705 

710 

715 

725 

745 

900 

100 mb 

.971 • 

.975 

.974 

.969 

.950 

.804 

.019 

-.620 

-.746 

-.801 

-.784 

-.748 

.970 

,974 

.973 

.967 

.944 

.762 

.181 

-.207 

-.345 

-.486 

-.488 

-.556 

150 mb 

.943 

.947 

.950 

.955 

.950 

.838 

.099 

-.562 

-.692 

-.745 

-.725 

-.686 

.944 

.948 

.952 

.954 

.946 

.800 

.243 . 

-.155 

-.293 

-.431 

-.432 

-.496 

200 mb 

.730 

.731 

.738 

.764 

.792 

.775 

.245 

-.348 

-.510 

-.579 

-.563 

-.533 

.736 

.737 

.746 

.769 

.792 

.730 

.272 

-.085 

-.226 

-.356 

-.357 

-.408 

250 mb 

-.528 

-.531 

-.524 

-.490 

-.433 

-.230 

.281 

.487 

.466 

.457 

.449 

.443 


-.534 

-.534 

-.528 

-.494 

-.441 

-.257 

.042 

.180 

.206 

.263 

.266 

.319 

300 mb 

-.870 

-.879 

-.879 

-.864 

-.839 

-.652 

.130 

.654 

.728 

.741 . 

.716 

.669 

-.868 

-.878 

-.878 

-.863 

-.838 

-.649 

-.118 

.206 

.311 

.422 

.419 

.475 

400 mb 

-.885 

-.895 

-.897 

-.889 

-.874 

-.708 

.084 

.650 

.742 

.763 / 

.738 

.683 

-.887 

-.897 

-.899 

-.892 

-.874 

-.694 

-.136 

.215 

.333 

.452 • 

.449 

.499 

500 mb 

-.839 

-.900 

-.902 

-.900 

-.887 

-.730 

.069 

.655 

.761 

.786 

.760 

.701 

-.894 

-.904 

-.907 

-.903 ■ 

-.887 

-.704 

-.129 

.235 

.364 

.486 

.482 

.529 

700 mb 

-.897 

-.90S 

-.906 

-.903 

-.877 

-.731 

.074 

.659 

.773 

.810 

.786 

.740" 

-.895- 

-.904 

-.904 

-.901 

-.881 

-.694 

-.109 

.255 

.386 

.518 

.514 

.571 

850 mb 

-.864 

-.870 

-.869 

-.868 

-.851 

-.712 

.061 

.625 

.750 

.803 

.784 

.754 

-.857 

-.864 

-.863 

-.862 

-.840 

-.668 

-.101 

.248 

.384 

.522 

•,522 

.592 


TEMPERATURE AND RADIANCE CORRELATION RATRIX FOR CLOUD FREE AND CLOUD CONTAMINATED DATA 
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Cloud Free Cloud Affected 


Tioo = -151.93 . 1.78 N^^^ 

Variance Explained = .950 

Standard Error of Estimate = 2.44 

Tioo = -151.53 . 1.78 N 323 . 

Variance Explained = .949 

Standard Error of Estimate = 2.46 

Ti50= -130.93. 1.54 N^^^ 
VE = .911 SEE =2.78 

Tiso= -131.20. 1.55 N^g^ | 

Ve = .911 SEE = 2.80 

T^Oo = -107.18 - 1.01 N ^^3 
VE = .675 SEE = 3.18 

T 200 = -91.27 - 0.92 N^25^1-54 N,g 3 
VE = .691 SEE = 3.08 

T 250 = -69.47 - 1.54 N^25 * 1‘06 
. 1.06 N ^33 . 0.97 N^oo 
VE = .506 SEE = 3.06 

T 250 = -48.95 - 0.72 N 323 - 0.90 N^g 3 

.1.26 N2.o-0.38 N 223 . 0.09 N... 
VE = .508 SEE =3.02 

T 30 O = \85 ^ '^705 

. 1.07 N 203 - 0.12 N 2,3 
VE = .835 SEE = 2.26 

■TSOO = - '-21 "685 " N 7 OO 

- "745 

VE = .802 SEE =2.54 

■^400 = "685 * 0-^8 N 205 

VE = .838 SEE = 2.68 

T 400 = 10.85 - 1.02 N 333 . 0.29 N 2.3 
VE = .820 SEE =2.85 

Ts^O = -7.32 - 0.85 N^g 3 . 0.60 N^^^ 
VE = .842 SEE = 2.62 

T 500 =20.81-1. 01 N^85. 0.39 N 705 
VE = .836 SEE =2.66 

T 700 = 5.82 - 0.42 N ^33 - 1.00 N^^^ 
.1.19 N^os 

VE = .862 SEE = 2.87 

T 7 OO = "685 ^ N 705 

VE = .837 SEE = 3.07 

Tgso = 33.98 - 0.94 N^^s ^ °.34 N^23 
VE = .786 SEE =4.12 

T85. = 61.11 - 1.16 N675 . 0.15 N 725 
VE = .765 SEE = 4.24 

1 


TABLE 6 

REGRESSION EQUATIONS FOR CLOUD FREE AND CLOUD AFFECTED DATA SETS 

Tin °C, N in erg/ 
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The most striking feature of Table 6 is the similarity between 
the cloud-free and cloud-affected regression parameters. There is 
hardly any difference in the variance explained between the two data 
sets and in some cases (200 and 250 mb) the variance explained is 
actually higher for the cloud contaminated set than for the cloud-free 
set. There is a significant difference, however, in the regression 
coefficients in the equation for 300 mb and below. In these equations, 
the coefficients of radiances such as N(705), N(710) etc, are larger 
in the cloud- free set than in the cloud-affected set. For example, 
for T^qq 'tke size of the coefficient for N(705) is greater in the cloud 
free set than in the cloud affected set. This difference is because of 
the effect of cloud on N(705). 

Both data sets contained data from 20-80°N. The greatest contribu- 
tion to the total variance of the temperature fields is the north-south 
variation. The variability of temperature around latitude circles is 
much smaller than the north-south variation. This leads to the high 
explanation of variance even for pressure levels in the low and mid- 
troposphere by equations containing radiances with weighting function 
peaks in the stratosphere (e.g., N(675) in equation for Tj,^_ ). This 
is because of the strong negative correlation between troposphere and 
stratospheric temperatures. Thus the regression equations in Table 6 
are probably of limited use, except perhaps for the top three levels. 

Below that the equations are explaining mostly the north-south temperature 
gradient and will not accurately pick out longitudinal, synoptic varia- 
tions in the troposphere. 
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2.3 Determinatioil of Tropopaiise Pressures and Tropopause Temperatures 

from IRIS Data 

2.31 Description of method 

In general the tropopause slopes from high altitude and low pres- 
sures over the tropics to lower elevations and higher pressures in 
temperate and polar latitudes. The temperature of the tropopause also 
varies with latitudes being cold in low latitudes and warmer in higher 
latitudes. Tropopause pressure and temperature also vary in relation 
to synoptic features with low and warm tropopause associated with upper 
troposphere troughs and cold and high tropopause associated with ridges . 
The IRIS radiances are related to atmospheric temperatures through the 
radiative transfer equation. IRIS channels with low wave numbers, for 
example with v = 667 cm ^ or v = 675 cm”^ , have a peak in their weighting 
functions (dx/^ Jlnp) in the low stratosphere while channels with higher 
wave numbers have weighting function peaks lower in the atmosphere. 

Thus, the radiances in the channels of the IRIS experiment are related 
to a weighted mean Planck function with different channels representing 
different overlapping layers. Because the Planck function is a function 
of temperature the weighted mean Planck function is related to the 
temperature in the layer. 

Figure 6 shows two plots of IRIS radiance against channel number 
with lower channel numbers associated with higher wave numbers. Channels 
with lower wave numbers which have weighting function peaks in the strato- 
sphere are higher on the ordinate than channels with weighting function 
peaks in the troposphere. The plots are for the same orbit on 3 July 1969 
for latitudes 24°N and 62°N and are for cloud-free data. The two pro- 
files are typical for low and high latitude radiance profiles. There is 
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Fig. 6 - IRIS Radiance Profiles 


t 


50 


a substantial difference between the two profiles. The low latitude 
profile (24.5°N) has higher radiances in chanels 1 through 4 and lower 
radiances in chanels 6 through 13. This cross-over of the two radiance 
profiles corresponds to a warm troposphere and cold stratosphere in 
low latitudes and a cold troposphere and warm stratosphere in high 
latitudes. The two radiance profiles also display a difference in the 
location of the point where there is a marked change in slope of the 
profile. On the low latitude profile the point of slope change is 
at channel 8; on the high latitude profile the slope change is at 
channel 6. The point at which the slope changes in the radiance profile 
is associated with the height of the tropopause. The tropopause is 
lower at 61.7°N than at 24.5°N and this is reflected in the lower 
position of the slope change at 24.5°N in Fig. 6. This example shows, 
at least in a general sense, that the position on the radiance profile of 
the change in slope is related to the tropopause position. The technique 
described in this section objectively determines the point in the radiance 
profile which is associated with the tropopause. 

The radiance profiles to be used in the tropopause computations, 
such as the profiles in Fig. 6, are determined by eliminating data 
points with unreasonable radiances, elimating data determined to be 
contaminated with clouds (using the 735 cm ^ channel) and averaging 
over five data points along the orbit. To the resulting radiance 
proriles a fifth degree curve is fitted by a least squares technique to 
channels 5 through 10 (705 cm ^ through 680 cm ^) to give an equation 
of the form 

N = a^ + a^^C + + . . . 


( 6 ) 
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where N is the radiance and C is the channel number which varies 
in this case from 5 to 10. 

To pick off the point in the radiance profile which is associated 
with the tropopause an indicator is needed of where in the profile the 
change in slope is greatest. After some subjective experimentation an 
indicator K defined by equation (7) was determined to do a good job 
in finding objectively the point the author selected subjectively. 

d^N 

K = dC^ (7) 



The tropopause channel , in other words the point along the C 
axis of the profile which is associated with the tropopause, is defined 
as the point where K is a maximum. The tropopause radiance (N.p) is 
defined as the radiance N in equation (6) when C = C.^.. 

A computer progrjim has been written to determine C.^, and N.^, 
from the IRIS data. The output of the program in addition to and 
are the position and time of the observation and the associated 
radiance profile. 

2.32 Con 5 )arison of satellite-derived parameters with tropopause pressures 
and temperatures obtained from radiosondes 

The parameters C.j, and obtained from the IRIS data should 

be in some way related to tropopause pressure and temperature respectively. 
The two parameters should represent features of a mean tropopause, i.e. 
a tropopause corresponding to a temperature profile smoothed with respect 
to height or pressure. The reason that the C.^, and obtained from 

the IRIS data should be related to a tropopause of a smoothed temperature 
profile is that the radiances themselves represent a smoothed structure. 



However, to objectively compare and with some actual 

observations of the tropopause, a set of near-simultaneous observations 
of and tropopause pressure and temperature were collected. 

The set of observations covers a 7 day period from 3 July to 9 July 
1969. The data set was confined to North America where tropopause 
pressure and temperature is given as a significant level report in 
the northern hemisphere data tabulations. The radiosonde observations 
used had to be within 180 nautical miles and six hours of the satellite 
position and time. Stations with double tropopauses were eliminated. 
Figure (7) shows a scatter diagram of tropopause pressure versus tropo- 
pause channel . There is a large scatter, but an apparent relation 
between the two variables with low associated with high tropopause 

pressure and high C.p with low tropopause pressure. The scatter is 
obviously not gaussian but skewed so that the long tail of the scatter 
is towards high pressure. This type of skewed scatter is evidenced 
by the points falling at higher pressure than the bulk of the points. 

These points represent reported tropopauses on the radiosonde data that 
are probably lower than the smoother profile tropopause that the parameter 
C.J, represents. A least square regression between tropopause pressure 
(P^) and is not applicable because of the non-gaussian scatter. For 

this reason the data set was divided into seven groups according to 
C.J, so that each group contained approximately the same number of data 
points. For each group a median tropopause pressure was calculated 
and plotted against the mean of the group. This plot is given in 

Fig- (8) with a curve fitted to the seven medians. From the fitted curve 
C.J, of 8.0 represents a tropopause pressure of about 115 mb, equal 

to 5.6 represents a of about 280 mb. 



5.4 5.6 5.8 6.0 62 64 6.6 68 70 7.2 7.4 76 7.8 8.0 8.2 

TROPOPAUSE CHANNEL Cj 


Scatter Diagram of Tropopause Pressure and Tropopause Channel 
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Fig. 8 - Group Medians of Tropopause Pressure versus Group Means of 
Tropopause Channel 
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Figure (9) shows the scatter diagram between tropopause temperature 
(T^) and tropopause radiance, N^. Again the data is separated into 
seven groups according to and medians in each group are calculated 

and plotted against mean tropopause radiance in Fig. (10). A linear 
fit to the plotted points is also given. An of 59 ergs/ is 

related to a tropopause temperature of -46°C, while an equal to 

42 ergs/ is about equal to a T.j, of -70°C. 

2.33 Comparison of analyses of tropopause channel and tropopause radiance 
with synoptic conditions at 200 mb 

Figures (11) and (12) show analyses of and N.^, respectively 

for 3 July 1969. The Cp analysis (Fig. 11) displays, in general, 
large over the tropics, low C.p in higher latitudes. On the 

analysis (Fig. 12), which corresponds directly to tropopause temperature 
(T.p) , low values are seen in low latitudes corresponding to cold tropo- 
pause temperatures, and high values in higher latitudes corresponding 

to warmer tropopause temperatures. Also evident from the diagrams are a 
number of synoptic features which appear on both charts. Low values 
of C.p which correspond to low tropopauses are related to high values 
of N.J, and therefore warm tropopause temperatures. The low and warm 
tropopause is a feature of troughs while high and relatively cold 
tropopauses are associated with ridges. In addition, tight horizontal 
gradients of tropopause pressure or tropopause temperature will be 
related to regions of a sharp sloping tropopause which in turn is 
associated with positions of jet streams. 

The troughs, ridges and jet streams that are associated with the 
features on Figures (11) and (12) can be seen in Fig. (13), an analysis 
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Fig, 11 - Northern Hemisphere Analysis of Tropopause Channel C,j, for 
3 July 1969 





Fig. 12 - Northern Hemisphere Analysis of Tropopause Radiance for 

3 July 1969. Units are ergs/. 
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Fig. 13 - Northern Hemisphere Analysis of 200 mb Geopotential Height Field 
from NMC Grid Data for 3 July 1969. Units are tens of meters. 
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of the 200 mb height field for 1200Z 3 July 1969, from NMC grid dta. 
Comparison of the analysis (Fig. 12) with the 2O0 mb analysis 

shows a good correspondence of features between the two charts. Troughs 
on the 200 mb chart coinicide with high radiance values associated 
with warm and low tropopauses. This association is evident, for example, 
in the trough centered about at 60°N, 70°W and its relatively weak 
extension eastward over southern Greenland and across the Scandinavian 
peninsula as seen in Fig. (13) . The trough and its eastward extension^ 
is associated with high in Fig. (12). 

The position of the jet streams is also marked on the anr-lysis 

(Fig. 12) . Over North America there is a confluence of lines 

with a relatively tight gradient eastward from about 110°W. The area 
of tight gradient narrows over the western Atlantic and extends eastward 
across the North Sea and the Paltic Sea in northern Europe. This 
position of the tight gradient coincides very well with the jet 

stream position as given by the 200 mb chart in Fig. (13). Tight 
gradients over the east coast of Asia and the western Pacific Ocean 
and over the eastern Mediterranean Sea and eastward are also associated 
with jet streams. 

Thus the tropopause channel C.^, , which is indicative of the tropo- 
pause pressure, and the tropopause radiance N.^, , which is associated 
with tropopause temperature show good agreement with synoptic features 
at approximate tropopause level. 

For the same day (3 July 1969) maps of C.p and N.^, were plotted 
and analyzed for the Southern Hemisphere with mixed results. In low 
latitudes, up to about 45°S, the analysis looks reasonable. At higher 
latitudes it becomes more difficult for the technique to determine the 



62 


tropopause information from the IRIS data. This difficulty is due 
to the presence of the stratospheric polar vortex in the higher latitudes. 
The cold stratospheric temperatures associated with the vortex reduce 

I 

the radiance in channels with a significant contribution from the 
stratosphere. This leads to a reduction of the curvature of the radiance 
profile and an increase in the difficulty of picking the point in the 
profile to represent the tropopause. Even with this drawback positions 
of troughs can be deduced and seem to compare well with satellite cloud 
observations. Further work will be done on this high-latitude, winter 
hemisphere problem in the next few months. 

Figure (14) gives the total ozone map for 3 July 1969 over the 
Northern Hemisphere. The features on the ozone map agree, in general, 
with both the features on the C.j, and maps in Figures 11 and 

12 and with the 200 mb chart in Fig. 13. This is to be expected in 
light of the work by Lovill (1972) including the determination of a 
relation between horizontal gradient of total ozone and jet stream 
speeds. Thus the gloval total ozone maps and the tropopause maps are 
comp 1 ement ary . 

2.4 Supporting Data 

In order to match and check results from the satellite experiments 
over the northern hemisphere, programs have been developed to read 
from tape and objectively analyze standard meteorological data archived 
each day by the National Meteorological Center, NMC. Special 
statistical characteristics of the atmospheric circulation (see the 
previous annual report) were compiled from rawinsonde data over both 
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hemispheres. These results, compiled in a M.S. thesis by Mr. R. Abbey, 
provide an indication of the relative importance of various planetary 
and cyclone scale wave regimes in the respective hemispheres. Thus, 
they aid and guide the satellite-based research. 

During the course of his research using ozone amounts inferred 
from IRIS, Dr. Lovill took full advantage of his extensive collection 
of ground-based ozone measurements. A library of satellite photo 
montages in the Atmospheric Science department is used extensively 
by the faculty and students for the Grant research. 

Near the close of this reporting period, the first HOLE data 
were received by the project. Tapes can be read and use of the data 
will be described in the next report. 
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3.0 PROGRAM FOR THE NEXT REPORTING PERIOD 

The ozone analysis pursued by Dr. Lovlll should be extended to 
additional seasons in both hemispheres. From the "ozone waves" 
which appear in a continuity diagram presentation we can establish 
phase-wind group velocity effects that can be used to arrive at 
planetary and cyclone wave lengths. The behavior of these waves 
during different seasons in the two hemispheres should give some 
indications of cyclone frequencies and intensities. Especially 
in the southern hemisphere these data can be tested against the 
findings from the EOLE balloon experiment. These comparisons can be 
made on a statistical basis only, because the period of the EOLE 
program does not fit the NIMBUS data period. 

The mapping of tropopause heights will be 'compared more closely 
with horizontal gradients in total ozone, both sets of data giving 
an indication of the position and intensity of jet stream. Attempts 
will also be made to use the information on tropopause location to 
improve the vertical temperature profiles derived from satellite data 
by inversion techniques. Instead of stratifying the data sample into 
latitude bands we will attempt to apply inversion methods to data 
arranged according to tropopause pressure and temperature as determined 
from IRIS data. 

Regression techniques, especially non-linear regression equations 
will be tested further as to their usefulness in arriving at the 
upper tropospheric and stratospheric temperature distributions from 
radiance data, and at geopotential height fields of constant pressure 
surfaces. Again, independent knowledge of the tropopause location. 
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derived from IRIS data might help to improve these regression techniques. 
The data base of geopotential height and temperature maps derived from 
satellite data will be broadened considerably beyond the data set, for 
which the regression method has been tested so far. 

With more analyses available, especially from the southern 
hemisphere, we will attempt to estimate the. energetics of the general 
circulation as outlined on p. 45 of our last Annual Report (Reiter, 
et al., 1971). Especially values of available potential energy and 
kinetic energy will be derived from satellite data and the significance 
of the variability with time will be tested. 

Our ozone analysis techniques, together with stratospheric 
temperature and geopotential height analyses derived from satellite 
data should be extended to NIMBUS 4 data. Results from these studies 
will lend themselves to comparisons with the Oxford analyses derived 
from the ^elective £hopper Radiometer flown on Nimbus 4. 

Analyses of the water vapor distribution, even though they were 
not pursued vigorously during the present grant period, will be considered 
as helpful augmentation in developing quantitative descriptions of 
flow patterns in the troposphere. 
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4.0 CONCLUSIONS AND RECOMMENDATIONS 

Following the early work of testing and program development, 
the NIMBUS-3 ozone analysis is now complete (see Appendix I) . Special 
areas for further study of the ozone data have been noted. Of course 
we also recommend that NIMBUS-3 results be checked using NIMBUS-4 
measurements. 

We have developed cloud checks independent of surface data, that 
will allow satisfactory study of atmosphere thermal structure in both 
hemispheres above the 700 mb level. We use inversion of radiances, 
definitely below 200 mb; regression on radiance for inferring the 
thermal structure at 100 mb and above. Calibration of our methods 
over the N. H. is satisfactory, once we have eliminated an apparent 
problem with the radiosonde temperatures reported at 30 mb. 

In addition, a promising and economical method to infer tropopause 
pressure and temperature from the IRIS radiances is available. It 
may aid the general thermal inversion and regression measurements. It 
definitely complements these analyses and also those of the ozone amounts. 

As this report is written we have already moved into the analysis 
and diagnostic phase of the thermal structure portion of the project. 
Development and testing has been slow but profitable. Meanwhile the 
data base for analysis (IRIS, SIRS) has been increasing. Fortunately, 
in addition to NMC data we will now have global SCR and southern hemi- 
sphere HOLE data to compare with our diagnostic results. 

We are most anxious to apply the developed techniques to a large 
data base to allow the most meaningful analysis of interhemi spheric 
circulation differences. If NASA can provide the basic data tapes, our 
experienced personnel (faculty, students and programmers) are ready 
for production runs on the CSU and NCAR computers. As in our prior 
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ABSTRACT 


Ozone is an important atmospheric trace constituent. The 
depletion of solar radiation between approximately 2000 and 3000 A 
is the result of strong absorption by ozone in the ultraviolet wave- 
lengths. The energy absorbed in this process is the prime source 
of thermal energy in the stratosphere. Because of this, ozone plays 
an important role in the large-scale motions of the atmosphere. As 
a result, there is an ever-increasing interest in global total ozone 
measurements. The paucity of total ozone measuring stations over 
vast regions of both the Northern and Southern Hemispheres has greatly 
hindered any complete understanding of the role stratospheric ozone 
plays in the general circulation processes of the atmosphere. With 
total ozone sensors on earth-orbiting satellites the ozone over the 
entire earth is remotely measured daily. The conclusions of some of 
the analyses from the first satellite (Nimbus III) are discussed in 
this paper. 

A strong correlation was found between the meridional gradient of 
total ozone and the wind velocity in jet stream systems. A procedure 
evolved which will allow the computation of the wind velocity at 250 mb. 
This should prove to be a useful relation to apply over the oceanic 
regions of the Northern Hemisphere and throughout most of the Southern 
Hemisphere. Wind velocities at jet aircraft operation levels could then 
be forecasted by consulting the daily satellite total ozone analyses. 
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A study of the total ozone distribution over two tropical storms 
indicated that each disturbance was associated with a distinct ozone 
minimum. 

A comparison of time-longitude stratospheric radiance values at 
60°S with values of the total ozone indicated that low (high) radiance 
values corresponded very closely with the low (high) ozone variations . 
The speed at which these ozone 'waves' progress eastward is greater 
in the winter hemisphere. The speed of eastward progression decreases 
as one approaches the lower latitudes in the winter hemisphere. In the 
equatorial region and in the Northern Hemisphere summer there is not a 
strong eastward progression of the ozone 'waves' but a westward 
progression. 

The average latitudinal distribution of ozone is described for the 
first time, both temporally and spatially, during the autumn and winter 
in the Southern Hemisphere. The lowest mean total ozone value for the 
globe is observed at 6°S. 

Maps are presented that describe the distribution of total ozone 
for the entire period (April to July 1969) for the Northern and Southern 
Hemispheres. Closely spaced total ozone isolines appear to reveal a 
rather strong jet stream west of South America at approximately 35°S. 

A large region of low total ozone was observed in the tropics, extend- 
ing from Southeast Asia eastward to the central Pacific. 

Both the Northern and Southern Hemispheres appear to provide 
'anchoring' mechanisms for total ozone ridges. The most predominant 
ridges are also associated with the 'tightest' ozone gradients. In 
the Northern Hemisphere these locations are primarily Eastern North 
America, Central Europe and Eastern Asia. In the Southern Hemisphere 
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ridges are evident over the East Indian Ocean, the Central Pacific 

Ocean and the Eastern Atlantic. The ridges correlate closely with 

upper atmosphere height contours. The ridge pattern is easier to 

delineate in the Southern Hemisphere winter than in the Northern 

Hemisphere summer. The 'tightest' gradient of total ozone in either 

hemisphere was seen over Japan. This correlates closely with the 

high frequency of jet streams observed over the region. 

An average global value of the total ozone for the entire period 

33c3 

was calculated to be SUB* m* atm* cm. 


James Edward Lovill 
Atmospheric Science Department 
Colorado State University 
Fort Collins, Colorado 80521 
February, 1972 
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1.0 Introduction 


1.1 The Nimbus III Spacecraft 

The Nimbus III spacecraft completes a sunsynchronous , circular 
orbit 1109 km above the earth once every 107.3 minutes. The satellite 
was launched into orbit at 07:54:03 GMT on 14 April 1969. Nimbus III 
represented a major upgrading of the Nimbus I and II spacecraft. 

Figure 1.1.1 represents the basic configuration of the Nimbus III 
spacecraft. The location of various experiments and spacecraft controls 
is indicated in the figure. 

One of the major experiment packages added to this satellite was 
the Michelson infrared interferometer spectrometer (IRIS) . The IRIS 
instrument has several scientific objectives. It has been providing 
information on the structure of the atmosphere and the emissive 
properties of the earth's surface. It is designed to measure the 
thermal emission spectrum of the earth between 5 and 25 microns 
(400 and 2000 cm with a spectral resolution of 5 cm (Spectra 
cannot be used between 5 to 7 p (1400 to 2000 cm because of de- 
creased responsivity of greater than thirty percent in this portion of 
the wave band.) In the spectral region covered, ozone, water vapor, 
and carbon dioxide bands are available in addition to spectral features 
associated with such minor atmospheric constituents as methane and 
nitrous oxide. The specific intensity measured in the absorption bands 
and in the transparent 'windows' can be used to determine vertical 
profiles of ozone, temperature, and water vapor. 

The essential nart of the interferometer (Hanel, 1969) is the beam 
splitter which divides the incoming radiation into two essentially equal 
components. After reflection on the fixed and moving mirrors. 
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Fig. 1.1.1. The Niinbus III Spacecraft. The locations of individual 
experiments are indicated. The Infrared Interferometer Spectrometer 
(IRIS) is at the lower right of satellite. 


respectively, the two beams interfere with each other with a phase 
difference proportional to the optical path difference between 
both beams. The recombined components are then focused onto the 
bolometer where the intensity is recorded as a function of path 
difference, 6. Siivce the mirror motion is phase locked to the 
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spacecraft clock, the mirror path difference Is also proportional to 
time. 


The central fringe may be light or dark depending on the 6 between 
the two beams. For polychromatic radiation, neglecting constant terms, 
the signal (1) at the detector, referred to as the interferogram, is 

i(6) = /k(B_ -B_ ) cos (itv6 -4>v) dv 1.1(1) 

The amplitude is nroportional to a responsivity factor, , and the 
difference in radiance between the scene within the field of view, B^ , 
and the instrument, B^. The phase is defined with respect to a point 
chosen as close as possible to, but not necessarily at, the zero path 
difference point. Imperfect optical compensation and residual phase 
shift in the analog part of the data spectrum cause the angle, 4>, to 
depend on the wave number. Reconstruction of the spectrum is by a 
ground based IBM-360 computer at the Goddard Space Flight Center (GSFC) . 

A sample spectrum calculated from a single interferogram is shown 
in Figure 1.1.2. One interferogram is obtained every sixteen seconds. 
Twelve and one-half percent of the IRIS measurement time is consumed in 
a calibration procedure. After 14 interferograms are taken in the oper- 


ating mode a comparison is performed with two known calibration sources 
to establish an absolute intensity scale. An onboard, warm blackbody 
provides one calibration interferogram and the interstellar void 
(assumed a nearly zero degree sink) the other (Nimbus Project, 1969). 

Interferograms are transformed by means of the Cooley-Tukey method 
(Cooley and Tukey , 1965). An anodisation function is applied to each 
interferogram prior to transformation in order to reduce side lobes 
of the instrument function. The particular apodisation function has 
the form (0.54 + 0.46 cos 2 ttt/T) where t is the distance from the 
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Fig. 1.1.2. The Tiifrareil 1 ii terferoineter Spectrometer thermal emission 
spectrum of the earth which was taken over the equatorial Pacific Ocean 
on 15 April 1969 (Nimbus Project, 1969). 


center and T is the total length of the interf erogram (see also, 
Blackman and Tukey, 1958). 

The Fourier transformation yields two amplitudes 

a = / ^ A(6) i (6) cos 2 itvddd 1.1(2) 

and 

b = A(6) i ((^) sin 2 nv6d6 1.1(3) 

where A(6) is the anodisation function. From a and b a magnitude 

V V 


C' 

V 



1.1(4) 
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and a phase 

<)) = arctan b /a 1.1(5) 

V V V 

are computed. The location of the pronounced central peak of the 
calibration interferograms is taken as the phase reference point . For 
atmospheric interferograms which may or may not have a strong central 
peak, the average phase reference point of several calibration 

interferograms is used (Hanel, 1969). This spectrum reduction is 
described in a more complete manner by Hanel (1969) , in which he 
delineates a method to eliminate noise components and defines the 
final calibration procedure. 

The field of view of the IRIS instrument is 4° half angle (Nimbus 
Project, 1969), This corresponds to an area on the earth's surface with 
a diameter of 150 km. 

High quality IRIS data were obtained until the experiment failed 
on 22 July 1969. Failure was attributed to the opening of the compen- 
sating clip at the bolometer which produced a large positive transient 
at the input to the bolometer pre-amplifier. One hundred days of data 
were obtained before this point in time, of which 95 days are used in 
this paper. 

1.2 The Photochemistry of Stratospheric Ozone 

The photochemistry of the stratosphere in general and of the ozone 
problem specifically was first thoroughly discussed by Chapman (1930) . 
Further development and refinement of the theory have been advanced by 
Diitsch (1946, 1956) and Craig (1950). 
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Many different approaches have been taken toward a better under- 
standing of upper atmospheric ozone photochemistry. Hunt (1966) has 
discussed the photochemistry of ozone in a moist atmosphere. Schiff 
(1964) and Norrish and Wayne (1965) have investigated the photochemistry 
of ozone in an oxygen-hydrogen atmosphere. Crutzen (1971) has dis- 
cussed the photochemistry of ozone in an oxygen-hydrogen-nitrogen oxide 
atmosphere. Complicated reactions like those discussed by Crutzen 
(1971) may number as many as fifty-six. 

To a first approximation it is sufficient to assume only an 
oxygen atmosphere for the discussion here. Ozone or triatomic oxygen, 
0 ^, may be formed by a three body collision 

0 + O 2 + M ^ 0^ + M. 1.2(1) 

M is any third atom or molecule. In order for the reaction in 1.2(1) 

S 

to occur, atomic oxygen must be formed. This will happen if a quantum 
of solar energy dissociates O 2 » 

0 2 + hv 0 + 0. 1.2(2) 

The reaction in 1.2(2) will occur only if A < 2423 A. 

Ozone may be destroyed in two ways: by collision, 

0 3 + 0 ^ 2 O 2 1.2(3) 

or by photo-dissociation, 

0^ + hA 0 + O 2 . 1.2(4) 

The reaction in 1.2(4) will transpire only when A < 11,000 A. 
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As one goes from the lower stratosphere to the mesosphere the 
ratio O 2 /O approaches zero. As a result there is less O 2 available for 
reaction 1.2(1) and ozone is not formed as rapidly. Additionally, there 
is more atomic oxygen, 0, available so that reaction 1.2(3) occurs more 
rapidly and thus ozone is destroyed. When these, and other reactions 
unfavorable to ozone formation, occur, the net result is a decrease of 
ozone. The height at which such a decrease usually commences in the 
middle latitudes is between 20 and 25 km. Although the ozone partial 
pressure may decrease above this level, the mixing ratio normally remains 
constant in the mesosphere at a value between 5 and 15 pg/g (see, e.g., 
Herlng, 1964; Lovill and Miller, 1968). 

Ozone in the lower stratosphere and upper troposphere is a quasi- 
conservative property and serves as an excellent tracer of atmospheric 
motions. This is clearly seen in the values presented by Paetzold (1953). 

He shows that the time required for ozone to decrease 1/e of its value 

2 3 4 

at 25 km is equal to 10 days; at 20 km, 10 days; and at 15 km, 10 days. 

Ordinarily the region near and at the surface of the earth will 
act as a sink for ozone. Little is known of the actual processes that 
bring about ozone decomposition in this region. Volz (1952) and Cadle 
and Ledord (1966) have shown that ozone can be decomposed upon contact 
with clouds, H 2 S, and other substances. On the other hand, investiga- 
tions by Haagen-Smit (1952) and Renzettl (1955) have shown that ozone 
also may be formed near the surface of the earth. This occurs when 
hydrocarbons in a polluted atmosphere are irradiated by sunlight . The 
direct result is large local concentrations of ozone. Therefore under 
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these conditions this region near the surface will act as a source 
rather than as a sink. 

The ozone concentration generally decreases in the vertical, 
from the surface to the vicinity of the tropopause, then increases to 
a maximum concentration at an average height in the middle latitudes 
of about 25 km, after which it steadily decreases with height. The 
height of maximum concenf-ation in the stratosphere varies considerably 
with latitude. The height varied from 14 km at 80°S (Lovill, 1965) 
during the Southern Hemisphere winter of 1963 to 26 km near the 
equator (Hering and Borden, 1964). 

1.3 Background of Early Investigations and Inadequacies 
of Present-Day Conventional Data 

Dobson (1930) was the first to initiate a large-scale total ozone 
sampling program. This began over a decade after Fabry and Buisson 
(1913) developed the technique for the basic measurements. For his 
complex of stations Dobson (1931) developed a spectrophotometer 
(usually referred to as the Dobson Total Ozone Spectrophotometer or 
simply The Dobson) for the measurement of total ozone. Over the past 
forty years the instrument has changed little and it is still essen- 
tially the same device originally used in Dobson's first network. 

With few exceptions the Dobson Spectrophotometer is the standard instru- 
ment for the measurement of total ozone in most countries. The major 
exception is the Soviet Union which uses the Filter Ozonometer (Bojkov, 
1969). 

Total ozone has also been measured by the integration of the ver- 
tical distribution of ozone determined from ozonesondes. The two 
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primary types of ozone sensors are the chemical (see e.g. Brewer and 
Milford, 1960; Komhyr, 1969; Lovill, 1969, 1970a) and chemiluminescent 
(Regener, 1960) devices. These ozonesondes are far more useful for 
their ability to measure vertical ozone distribution than for their 
total ozone measurements. This fact will become abundantly obvious in 
Chapter 3 when vertical and horizontal transport processes in the 
vicinity of a jet stream system are discussed in relation to early 
work by Briggs and Roach (1963) . 

Climatological studies of total ozone utilizing the Dobson Spectro 
photometer and the ozone sensors reveal the following for the Northern 
Hemisphere (for a summary see Ddtsch (1969)) . Total ozone reaches a 
maximum in the late winter and early spring. Minimum values occur 
during late summer and early autumn. The highest total ozone concentre 
tlons are recorded in the higher latitudes and the lowest values at the 
lower latitudes. Considerably less is known of the total ozone distri- 
bution over the Southern Hemisphere due to a serious lack of observing 
stations. Total ozone variations in the mid-latitudes of the Southern 
Hemisphere appear to be similar to those in the Northern Hemisphere, at 
least when data from sites in Australia are analyzed (most of the 
Southern Hemisphere large-scale ozone analyses have been conducted by 
Kulkarni (1966) and Pittock (1968)). 

The paucity of total ozone measuring stations over the oceans and 
in the tropics in the Northern Hemisphere and throughout the Southern 
Hemisphere has greatly hindered any complete understanding of the role 
stratospheric ozone plays in the general circulation of the atmosphere. 
Considerable information is available on the latitudinal variation 
of total ozone over the continents in the Northern Hemisphere (see 
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London, 1963) and the temporal variations of the vertical distribution 
of ozone at individual stations (Bering and Borden, 1964; 1965a, b; 
1967). But essentially no information is available over the Northern 
Hemisphere oceans. Little can be advanced about the longitudinal 
variations of total ozone except over the continents in the Northern 
Hemisphere. In addition, due to the sparsity of stations in the 
Southern Hemisphere, almost no understanding is available of the inter- 
hemispherical exchange processes dealing with ozone. 

The Nimbus III sensor provides two or three thousand total ozone 
data points over the entire globe each day. The analysis of these data 
will provide significant advances in research related to both ozone and 
the general circulation. These data will be discussed in the following 


sections . 
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2 . 0 Error Analysis and Data Preparation 

Data obtained from an earth-oriented satellite must be calibrated 
in some manner. This is usually accomplished by comparison with 
ground- truth data. In Chapter 2.1 the pertinent data readout 
obtained during a typical orbit is briefly discussed. Chapter 2.2 
discusses the merits of satellite data processed by the inversion 
method and the regression technique. In Chapter 2.3 a comparison is 
made between various surface ozone station data and data obtained by 
the Nimbus III satellite. A brief description of the method by which 
the data are computer-processed is given in the final section. 

2.1 Typical Data Orbit 

Theoretically, approximately 4000 atmospheric spectra could be 
obtained each day from the IRIS unit if all spectra were useful and all 
orbital data were transmitted to the earth receiving stations. In 
reality, however, only about 50% of the spectra have been reduced via 
regression coefficients to total ozone data. Noise problems within 
individual spectra eliminate many spectral data; additionally, not all 
data stored onboard the satellite are transmitted (see Chapter 4) . 

Dr. C. Prabhakara (GSFC, NASA) provided the computer tape with the 
regression coefficient ozone data on which the analyses in this paper 
are based. 

Figure 2.1 is representative of the typical total ozone distribution 
observed during a Nimbus III orbital ’sweep’ from the high latitudes 
of the Northern Hemisphere to the high latitudes of the Southern 
Hemisphere. The data shown in Figure 2.1 were obtained during orbit 
1077 on 3 July 1969. During this satellite pass the data orbit ’swept’ 
from western Siberia southward across the western Pacific and into the 




Figure 2.1 Total ozone distribution during orbital pass 
Number 1077 on 3 July 1969. Data have been 
analyzed by the Regression Technique. 


Southern Hemisphere across central Australia and ended near the coast of 
Antarctica after covering the eastern Indian Ocean. 

Even a cursory examination of Figure 2.1 indicates a large 
range in total ozone concentration from the polar regions to the 
equatorial area. Total ozone values approximately double from the 230 
m. atm. cm (milli* atmospheric -centimeter) amount near the equator to 
460 m-atm-cm at 58 S. In chapter 1 total ozone amounts on the average 
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were observed to decrease smoothly from high to low latitudes . This 
obviously is not the situation during this orbit. In fact the total 
ozone amount never decreases smoothly from the polar regions toward 
the equator. There are always perturbations among the best line fit. 

These perturbations as seen in Figure 2.1 are mostly the result of the 
redistribution of total ozone in the atmosphere by cyclones and anti- 
cyclones. The largest perturbations along the line in Figure 2.1 are 
associated with the strongest anticyclone/cyclone systems. Since over 
ninety percent of the ozone in the atmosphere is found in the stratosphere, 
it becomes quite obvious that the redistribution process will be most 
efficient when these systems extend to their greatest vertical extent 
into the stratosphere. If this reasoning is carried one additional 
step farther, then one may stipulate that when large height differences 
exist between antlcyclonic and cyclonic systems in the stratosphere, 
large changes in ozone will occur along this height gradient. 

This argument can be validated when the ozone perturbations in 
Figure 2,1 are compared with upper-atmosphere synoptic features. In 
Figure 2.1 the largest gradient of total ozone occurs near 30°S. At 
this location the satellite is measuring total ozone over south 
Australia. A strong anticyclone is evident over north central Australia 
and a cyclone is seen just off the south vest tip of Australia, both 
extending into the stratosphere. The height gradient between these two 
systems at 300 mb is the greatest found on earth on this particular date. 

A 'tight' gradient such as seen in this case between two pressure systems 
is indicative of a strong geostrophic wind velocity. This suggests 
a strong positive correlation between the geostrophic wind and the 
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total ozone gradient. This correlation will be explored further in 
chapter 3. 

Two other features, in Figure 2.1, both pertaining to the geogra- 
phical and seasonal bias of ozone distribution, are worthy of mention. 

The winter hemisphere as a greater amount of ozone in its atmosphere 
than the summer hemisphere (see chapter 1). Since Figure 2.1 represents 
the v/inter season in the Southern Hemisphere, one expects to observe 
more ozone in the Southern than in the Northern Hemisphere. This is 
indeed the case. There is 19 percent more ozone at 70°S than at the 
comparable northern latitude; the largest difference indicates 11 
percent more at 60°S than at 60°N. A second feature of this particular 
orbit indicates the global minimum of total ozone to be located in the 
Northern Hemisphere (the summer hemisphere) at 8°N. 

2.2 Comparison of Inversion Data with Regression Data at Daily Intervals 

Prabhakara (1969) , Prabhakara and Conrath (1970), and Prabhakara, 
et al. (1971 b,c) have discussed the technique by which total ozone 
concentration may be determined from the IRIS instrument on the Nimbus 
III satellite. Prabhakara and Conrath (1970) have determined a method 
to invert the radiance measurements made by the IRIS. A number of 
calculations were made using this technique to determine total ozone. 

The inversion results yielded a standard error of estimate of six 
percent when compared with the surface Dobson spectrophotometer data. 

A complete discussion of the ozone inversion technique and the 
associated errors can also be found in Russell (1970) . 

The inversion technique is a time-consuming computer method and 
therefore a costly means of determining total ozone on a global scale. 

For this reason regression coefficients were determined by Prabhakara. 
et al. (1971a) by which total ozone amounts could be produced. 
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These coefficients were based upon data obtained from April 1969 
Nimbus III readout. Table 2.2 is a comparison of the regression and 
the inversion methods for determining total ozone at 9.6p. All global 
total ozone data were compared for 23 April and 3 July 1969 . Since the 
April data were the base for determining the regression coefficients, 
the real test of the coefficients arrives when the comparison is made 
of the two methods on 3 July. However, the two dates indicate (see 
Table 2.2) an identical standard error of estimate of 3.7 percent when 
all latitudes are compared. As expected, in general, the error between 
the two techniques increases as one goes from low to high latitudes. 

This is because the remote sounding techniques do not perform well 
where vertical temperature profiles are nearly isothermal, as they 
tend to be in the lower stratosphere of high latitudes. 

Fifty-eight percent of the inversion data points were greater than 
those found from the regression coefficients for the 23 April comparison. 
The 3 July case, however, indicated that 53 percent of the regression 
data values were greater than the inversion technique values. It, 
therefore, appears as though there is little bias toward positive or 
negative values . 

2.3 A Comparison of Dobson Data with Regression Data at Monthly Intervals 

The Dobson Spectrophotometer measurements of total ozone are usually 
considered to be highly accurate. On an absolute scale the Dobson 
instrument probably has an absolute error no greater than one to three 
percent. In the following paragraphs these data are compared with the 


satellite data. 
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Prabhakara (1969) and Prabhakara and Conrath (1970) expended 
considerable effort to determine the regression coefficients to be used 
to convert the IRIS spectra to meaningful ozone data. Sufficient 
satellite total ozone data are available such that a study may be made 
which compares the average Dobson data with the average regression 
determined ozone values for the period . 

In order to have a comparison of the highest quality » it was 
necessary to eliminate some of the days during which a less than 
sufficient number of ozone data points were processed, from the 
satellite. Seventy-five days were selected and a map of the average 
global ozone distribution during these 75 days was constructed. 

Only 38 surface stations around the world took total ozone measurements 
during a sufficient number of days for significant statistical 
comparison with the IRIS data. Only six of these stations were in the 
Southern Hemisphere. Table 2.3 describes these results. 

Table 2.3 

Standard error of estimate between the average 
Dobson Spectrophotometer total ozone and IRIS 
regression coefficient determined total ozone 
values at selected locations. 


Number of days 
Dobson obs taken 

Number of sta- 
tions in sample 

Standard 

deviation 

74 - 75 

9 

3.6% 

73 - 75 

11 

3.9 

70 - 75 

19 

4.6 

65 - 75 

35 

6.0 

60 - 75 

38 

6.1 


Table 2.3 Indicates that when data were available on every day, or 
every day except one, that the two data sets compare quite favorably. 
The standard error of estimate of 3.6% is well within the acceptability 


limits . 
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When all 38 stations are included in the data sample this increased to 
6.1%. However, with this number of stations, as much as 20% of the 
days may be missing in some cases. 

The largest standard errors were found over the Soviet Union. 

The standard error in several instances exceeded 15% . In all cases 
the Soviet Filter Ozonometer reported a lower value of ozone than the 
IRIS instrument (see also Bojkov, 1969). 
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3.0 Total Ozone Gradients Near Strong Barocllnic Zones 

Reiter (1971) suimnarized past research work on the ozone distribu- 
tion in the atmosphere that could be derived from ozone soundings, air- 
craft measurements and ground-based umkehr techniques (For specific 
details see Lovill, 1965, 1969, 1970a, b, c; Lovill and Miller, 1968; 
Bojkov and Lovill, 1969; Reinking and Lovill, 1971). From these studies 
it becomes quite obvious that the dynamically induced vertical motions 
associated with jet maxima near the tropopause extend into the lower 
stratosphere and Induce an appreciable vertical transport of ozone. 
Specifically, it was found that the stratospheric air on the cyclonic 
side of the jet stream, which is subject to sinking motions, abounds in 
atmospheric ozone; whereas the warm air in the jet core and on the anti- 
cyclonic side of the jet maxima, being of tropospheric origin, is defi- 
cient in ozone. Since the ozone surplus on the cyclonic side, as well 
as the deficit on the anticyclonlc side, extends over a considerable 
depth of the atmosphere, one would expect to find a marked horizontal 
gradient of total ozone to be present across the axes of well-developed 
jet streams. 

Thirty-three cases of relatively sharp meridional gradients of 
total ozone in the Northern Hemisphere on the objectively analyzed charts 
of ozone as observed by Nimbus III, have been tested with respect to 
the presence of absence or jet maxima. Not only was a strong correlation 
found between such gradients of ozone and the proximity of jet streams, 
but the magnitude of these gradients appeared to stand in direct 
relationship with the intensity of the jet stream, i.e., with the peak 
wind speeds observed in each case. 
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In one individual case of 7 May 1969, a strong Pacific cyclone 
and its associated jet are shown superimposed on an objective analysis 
of the IRIS ozone values. Maxima and minima of ozone correspond well 
with the postulated vertical motion fields near the jet core. 

A tropical storm in the eastern Pacific and a hurricane in the 
western Pacific are shown with the objective ozone analysis superimposed. 
Both cyclones are associated with ozone minima. Horizontal and vertical 
circulation mechanisms of each are discussed in relation to total ozone 
distribution. 

3.1 Meridional Gradients of Ozone Near Jet Maxima 

A relation between total ozone gradients and jet stream wind maxima 
was first suspected when the daily maps of total ozone from Nimbus III 
were compared with the National Meteorological Center (NMC) wind 
analyses. Further work indicated that the strongest total ozone gradient 
was directly related to the highest wind speeds from the NMC analyses. 
Figure 3.1 represents the slightly non-linear relation that was found 
to exist between these two parameters after 33 case studies were 
plotted . 

The plot of the total ozone gradient vs. the jet stream core 
velocity was found to have the least scatter when winds at the 250 mb 
level were used. The plot represents IRIS total ozone and NMC wind 
velocity comparisons for areas throughout the Northern Hemisphere, from 
20°N to 70°N. The majority of the data were obtained between 30° and 
60°N, over the large land masses of North America, Europe, the Soviet 
Union, and China as well as the oceanic regimes of the north Pacific and 


Atlantic . 
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Fig. 3.1. Graph of the non-linear relation between total ozone 
gradient (in m. atm •cm/ km) and the jet stream core velocity 
(in ms~^) . (Numbers identify individual cases of study.) 
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There is little doubt, after even a brief study of Figure 3.1, that 
the dynamics of the jet stream are extremely effective in redistributing 
the total ozone in the stratosphere. It is generally felt (Reiter, 1963; 
Lovill and Miller, 1968) that on the cyclonic side of a jet stream 
system where air is descending, the stratospheric air is richer in 
ozone concentration than on the anticyclonic side of the jet maximum 
where ascending motions are prevalent and air poor in ozone is found. 

3 . 2 The Total Ozone Distribution Over a Strong Pacific Cyclone and the 

Associated Jet Stream System 

Figure 3.2 indicates in a remarkable manner the dynamics of a 
strong jet stream system, such as discussed earlier. In this figure 
(7 May 1969) a strong Pacific cyclone and its associated jet are 
shown superimposed on an objective computer analysis of the IRIS total 
ozone values. Several features are worthy of mention. 

As previously discussed, the stratospheric air on the north 
(cyclonic) side of the jet stream is much higher in ozone than the 
air on the south (anitcyclonic) side. In Figure 3.2 the large ozone 
maximum of 450 m*atm*cm to the northwest of the jet axis and of the 
cloud band is produced by the vertical motion around the jet axis and by 
horizontal advection of higher ozone amounts from the north. At the 
same time the circulation around the cyclonic system is responsible for 
the large ozone minimum centered at 40°N and 170°E. This minimum is a 
result of ascending air to the south of the system as well as the 
advection of ozone 'poor' air from the tropics. Notice the extremely 
strong north-south gradient of ozone isolines along longitude 120°E. 

From 27°N to 53°N there is an increase in ozone by 180 m’ atm -cm. 
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This is an increase in total ozone of over 60% in 26 degrees of lati- 
tude. The highest wind velocity along the jet is associated with the 
very closely spaced ozone isolines at 45 to 50°N near 160°E. This 
corresponds correctly with the findings in Section one. 


40 ms ' ISOTACH 

CLOUD BAND 

JET AXIS C = CYCLONIC CENTER 



Fig. 3.2. Total ozone contours (m.atm.cm) over a large Pacific 
Ocean extra-tropical cyclone and its associated jet stream (7 Mav 
1969) 
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3.3 Tropical Storm Ava and the Meso-Scale Total Ozone Distribution 
Tropical Storm Ava was first observed as a closed circulation 
by ESSA 9 off the coast of Guatemala on 1 July 1969. The location 
of the storm for the next several days is indicated in Figure 3.3. 



Fig. 3.3. Total ozone analysis (m- atm- cm) over Tropical Storm 
Ava on 7 July 1969. Ava indicated by A and the numbers 3-6 
indicate locations earlier in July. 


Ava reached a maximura intensity of 29 ms ^ just before a sharp turn to 
the north began on 5 July (see Figure 3.3). Ava had decreased in 
intensity by 7 July. The ESSA 9 cloud picture of the area for 7 July 




-as- 


indicated an inner closed cloud band of approximately 120 km diameter 
and an outer closed band of several hundred kilometers extent. The 
outer band is indicated by dashed lines in Figure 3.3. The storm 
appeared to have clouds associated with its movement that enclosed an 
area of 7° latitude by 7° longitude. 

The isolines of total ozone measured by Nimbus III during this j 
period have been superimposed upon the outline of this tropical 
storm which is in its dissipation stage. Penn (1965, 1966) has provided 
information on the distribution of ozone at various levels in two 
hurricanes. U-2 aircraft observations over one of the hurricanes 
indicated that from about 120 to 105 mb the region over the eye was 
warmer and considerably richer in ozone than other areas in the 
vicinity of the storm. Higher in the stratosphere (100 to 50 mb) 
the horizontal gradients were weak. 

The major difference between Penn's data and the data presented 
here is that whereas his measurements were of the local ozone dis- 
tribution at specific locations. Figure 3.3 represents the total ozone 
distribution over a large region which surrounds the tropical distur- 
bance. (A similar analysis of the total ozone distribution over 
Hurricane Bernice was undertaken by Allison, et al. (1971) using 
Nimbus III data.) 

In Figure 3.3 a broad region of low total ozone concentration is 
seen surrounding tropical storm Ava. The center of low total ozone 
is just to the west of the storm. The distribution of total ozone 
around the tropical disturbance suggests, weak though it may be, 
that there is sufficient rising motion and advection of lower latitude 
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air in the vicinity of the storm to account for the total ozone 
minimum. 

Additional investigations of the distribution of total ozone 
around tropical disturbances should shed more light on the stratospheric 
motion patterns associated with these storms. A stronger tropical 
disturbance will be discussed in the next section and it will be seen 
there that the total ozone is also a minimum when the storm is well 
above minimum hurricane strength. 

3.4 Typhoon Susan and the Meso-Scale Total Ozone Distribution 

Typhoon Susan began forming two days after the Nimbus III space- 
craft was placed into orbit. The total ozone data in Figure 3.4 were 
obtained on 22 April 1969, the third day after the IRIS unit became 
fully operational. Although the disturbance was first noticed on 16 
April, it did not attain typhoon strength until 20 April, 

On 22 April winds near the eye were in excess of 50 ms At this 
time the ESSA 9 cloud photographs indicated that the typhoon covered 
a very large region of the Pacific. Cloud bands extended from 3 to 
18°N and from 122 to 140°E, an area of fifteen by eighteen degrees. 

The two most obvious cloud bands visible from ESSA 9 are shown 
in Figure 3.4. The location of Susan on 22 April is indicated by the 
S in the figure. The movement of the storm during the past four 
days is also indicated. The typhoon has slowed its forward movement 
and by the 22nd has begun curving to the northwest. At the time of 




-27- 


its location in Figure 3.4, Susan is only 750 km southeast of Manila 
and has just begun to move onshore over the Philippine Islands. 



Fig. 3.4. Total ozone analysis (m. atm* cm) over Typhoon Susan 
on 22 April 1969. Susan is indicated by S. The numbers 18 - 21 
Indicate locations earlier in April. 


Figure 3.4 indicates that a distinct total ozone minimum is 
surrounding the typhoon. The center of the minimum is approximately 
390 km southeast of the typhoon center. This was the approximate 
location of Susan 48 hours previously. 
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It is suggested that the same mechanisms that were called upon to 
produce the total ozone distribution in tropical storm Ava may also 
be responsible for the total ozone distribution around Susan. The 
basic oddity of these two tropical storms is the following. i.he 
weaker tropical disturbance, Ava, had just as strong a total ozone 
minimum associated with it as did the stronger typhoon. One assumes 
that the vertical structure of the typhoon extends to considerably 
greater heights than that of a dissipating tropical storm such as 
Ava. Therefore the redistribution of total ozone around the more 
intense storm should be more obvious. This is not true in the case of 
the two tropical cyclones just examined in Figures 3.3 and 3.4. 
Evidently a great deal more study needs to evolve around this problem. 
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4.0 Global Ozone Climatology; Temporal and Geographical Variations 

4.1 General Remarks 

Dynamically induced subsiding motions, in combination with strong 
horizontal advection, behind a short-wave trough could produce large 
positive changes in the total ozone concentration over an extensive area. 

As was illustrated in Figure 3.2, the ozone concentration decreases 
rather markedly on the west side of a ridge: With these correlations we 

can interpret longer term temporal and geographical variations of 
ozone on a quasi-seasonal basis. Many general circulation character- 
istics may thus be portrayed by the ozone distribution. 

4 . 2 Northern and Southern Hemisphere Comparisons 

Several studies in recent years have shown distinct differences of 
the atmospheric general circulation of the Northern and Southern Hemi- 
spheres (for a summary, see Reiter, 1969). Whereas, earlier publications 
postulated a more zonal flow in the Southeim Hemisphere supposedly be- 
cause of less land mass and fewer mountain massifs, more recent investi- 
gations show that 'such statements may have to be qualified somewhat. 

Wooldridge and Reiter (1970) have shown that significantly stronger 
horizontal anisotropy of flow prevails at cyclone wavelengths in the 
Southern than In the Northern Hemisphere, with meridional perturbations 
exceeding zonal flow perturbations. Eulerian spectral densities of the 
zonal component in the Southern Hemisphere were half of those in the 
Northern Hemisphere. This would imply that there is less energy at 
cyclone-wave frequencies, u, in the Southern than in the Northern Hemisphere. 
Wooldridge and Reiter (1970) suggest that this occurs because of a 
relative lack of orographically induced planetary long waves in the 
Southern Hemisphere that are superimposed upon the hemispheric jet 
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stream zones. They concluded that the relatively zonal character of the 
Southern Hemisphere pertains to the long planetary-scale waves only. 

Starr and Dickinson (1963) have shown that standing eddies are very 
important in transport processes of the Northern Hemisphere. Obasi (1963a) 
pointed out that transient eddy processes are of utmost importance in the 
flux of momentum toward the pole in the Southern Hemisphere. Obasi (1963b) 
also found that the conversion of eddy kinetic to zonal kinetic energy 
by such waves was about twice as great in the Southern Hemisphere as in 
the Northern Hemisphere. 

These studies imply that transport processes in the Southern Hemi- 
sphere depend more on the cyclone scale disturbances than in the Northern 
Hemisphere. In the Northern Hemisphere much transport is associated 
with standing eddies which have wave numbers of 1, 2, and 3. 

Jet stream patterns in the Northern Hemisphere are strongly influ- 
enced by the hemispheric wave number three (Reiter, 1963). This results 
in the formation of the three strongest mean jet maxima over Japan, 
eastern North America, and Europe, (Kao and Hurley, 1962). On the 
average these jet stream systems are stronger in the Northern Hemi- 
sphere winter than in the summer. 

In the Southern Hemisphere Van Loon (1964, 1965) finds a differ- 
ently structured jet stream system: In winter, the zonal momentum is 

greatest. The momentum is distributed over two hemispheric jet stream 
systems. In the Southern Hemisphere summer only one jet stream system 
is reported (over the Indian Ocean) by Van Loon (1964) . This single 
jet maximum has a greater velocity than either of the two jet maxima 
observed during the winter in the Southern Hemisphere . 


Large-scale circulation patterns also show a disparity between the 
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Northern and Southern Hemispheres. Planetary waves number two and three 
shape the polar vortex in the Northern Hemisphere winter. In the 
stratosphere , the circulation is such that upward motion is seen over 
the North Polar Regions during the late winter during and after the 
breakdowns of the stratospheric polar vortex CReed et al, 1963; Teweles, 
1964) . Subsiding motions are wide-spread over the middle latitudes 
(Reiter, 1969) . Since ascending currents dominate in the equatorial 
regions, the general circulation of the stratosphere is accomplished by 
two cells which produce an observed spring maximum of ozone in middle 
and high latitudes of the Northern Hemisphere (Reiter, 1971) . London 
(1963) in a state-of-the art paper presented seasonal and annual mean 
charts of total ozone data from sixty stations, most of which were in 
the Northern Hemisphere. He found in both the seasonal and annual means 
that three ozone maxima were evident in troughs extending from the high 
latitudes toward the equatorial regions. These waves in the ozone iso- 
lines conformed quite closely to the troughs in the stratospheric cir- 
culation pattern. 

4.21 Inter-hemispherical Comparisons of Monthly Ozone Variability 

Figures 4.211.1a, 4.212.1b, 4.212.1c, and 4.212. Id represent computer 
analyses of Northern Hemispheric total ozone data and Figures 4.213.1a 
to 4, 213. Id show similar analyses for the Southern Hemisphere. The 
charts for April represent the last twelve days of the month and the 
analyses for July, the first twenty-one days of the month. More data 
were available for analysis in the Southern than in the Northern Hemisphere. 
For example in May in the Northern Hemisphere 22,519 points were analyzed. 

In the Southern Hemisphere 30,110 data points were available - an increase 
of 25%. This deficiency of points in the Northern Hemisphere is mostly 



Fig. 4.211.1a Total ozone distribution in the Northern Hemisphere 
for April 1969. (Total ozone values in m- atm* cm) 

due to data retrieval from the satellite at the two acquisition 
stations. All orbits cannot be recovered owing to memory storage on 
the satellite and to an insufficiently rapid data bit transmission 
from the spacecraft to the station. 
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4.211 A Detailed Study of April In the Northern Hemisphere 

Figure 4.211.1a shows the mean total ozone distribution for April 
in the Northern Hemisphere. Two features are obvious in this figure: 

(1) a 'tight' ozone gradient extends from eastern North America to 
about 60°W in the western part of the Atlantic Ocean, and over eastern 
Asia and Japan; (2) there are four to five major waves in the ozone 
isolines. 

An analysis of the National Meteorological Center (NMC) upper 
troposphere and lower stratosphere wind data indicates that the two 
regions with closely spaced ozone isolines coincide with the average 
location of strong baroclinic zones. The 'tight' gradient of mean 
total ozone concentration seems to be directly identifiable with the 
mean position of strong polar jet stream systems. 

Four or five long waves also can be Identified in the ozone iso- 
line pattern. These ozone ridges (maxima) and troughs (minima) cor- 
respond closely to long-wave features found in the height contour fields 
at 300 to 100 mb. Shorter, transient total ozone waves tend to be fil- 
tered out in the averaging processes used to obtain monthly means. 

Several of these short-wave features are evident in this "April" map 
(Fig. 4.211.1a) however, since only twelve days of data were averaged. 

The other months have two to three times as many days to average, and 
thus identify mostly standing wave features. 

In Figure 4.211.1a the areas from the 300 m* atm* cm total ozone 
isoline southward (northward in the Southern Hemisphere; see corresponding 
diagrams) of the thirty degree latitude line are widely dotted to in- 
dicate preferred global regions of southward total ozone transport by 
the standing, long planetary waves. The densely dotted areas from the 
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300 m- atm* cm isoline northward (southward in the Southern Hemisphere) 
of the thirty degree latitude line represent preferred global regions 
of northward total ozone transport. A large region over which the 300 
m- atm* cm line is displaced north of 30°N extends from Central Asia to 
western Japan. With the exception of slight displacements to the north 
over Iran and Arabia and the Mid-Atlantic, the remaining area of the 
Northern Hemisphere has the 300 m* atm. cm isoline displaced to the south 
of 30°N. This implies that a larger amount of ozone is advected south- 
ward of 30°N than is advected northward at this time of the year. The 
orientation of the 300 m. atm* cm contour is mostly a function of the 
standing planetary wave system and of the resulting subsiding (ascending) 
and southward (northward) mass transport which results in preferred 
regions of total ozone maxima (minima) . The average concentration of 
ozone along a latitude circle, however, is mostly a function of the 
zenith angle of the sun. 

The 260 m* atm* cm contour for the months of April-June, and the 
280 m* atm* cm contour for July have been marked especially, to emphasize 
preferred regions where ozone maxima exist in the tropics. In April 
(Fig. 4.211.1a) the 260 m- atm* cm contour dips far south and generally 
intercepts the equator over much of North Africa. 

4.212 A Comparison of General Features in the April to July Mean Maps 

In May, June, and July (Figs. 4.212.1b, c, d) a very ’tight' ozone 
gradient is evident again over eastern North America. It is worthy to 
note that this area of concentrated ozone contours has moved a consid- 
erable distance northward from April to July. In April the average value 
of total ozone off the east coast of North America at 45°N was 400 


m* atm. cm. By July the value at this location had decreased to 340 
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Fig. 4.212.1b Same as 4.211.1a except for May 1969. 

m-atm*cm. This seems to correlate reasonably well with the weakening and 
movement toward the north of the baroclinic ozones. The tightly packed 
ozone contours seen over Central Asia and Japan in April appear to per- 
sist through July, but with considerable weakening after June over Japan. 
By comparison the gradient over eastern North America in July is 'tighter' 
than that of the gradient over Japan during that month, but about equal 
to that found over eastern Asia. Even though there are four standing 



- 36 - 


100 



80 

Fig. 4.212.1c Same as 4.211.1a except for June 1969. 

planetary waves indicated in the height contour field (Scherhag, 1969c) 
and in the ozone contour field in July (Fig. 4. 212. Id), it appears as 
though much of the jet stream activity is concentrated over eastern 
North America and eastern Asia. The ozone gradients over the rest of 
the Northern Hemisphere are considerably less than over these two regions. 

Scherhag (1969a) reports that the average wind velocity for several 
Japanese stations between 35 to 40°N at longitude 135°E for May at 200 mb 
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Fig. 4. 212. Id Same as 4.211.1a except for July 1969. 

was 40 ms The total ozone gradient in this region, using Figure 
4.212.1b, was 0.54 in m.atm*cm/km. Using Figure 3.12.1, correlating 
ozone gradient and jet stream winds, we arrive at an estimated velocity 
of 37 ms This is in good agreement with the actually measured winds, 
even though the wind reports pertain to 200 mb, and the correlation chart 
was developed for use at 250 mb. 

From April to July in the Northern Hemisphere (Figs. 4.212.1b~d) 
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there are three to five long waves visible in the total ozone contours. 

The month of June is the only month with a dominant three wave pattern. 

A closed center of total ozone appears in high latitudes at approximately 
175°E longitude. A dominating feature on the 200 mb mean June contour 
map is a closed warm center west of Alaska (Scherhag, 1969b) in the same 
location as the closed ozone contours. This feature remained quasi- 
stationary during June and appears to have influenced rather markedly 
the upper atmospheric flow over western North America and northeast Asia. 

In April we find most of the 300 m*atm-cm contour to the south of 
latitude 30°N. By July (Fig. 4. 212. Id) much of this contour lies to the 
north of this latitude. The months in-between represent a transition 
period during which the general circulation in the Northern Hemisphere 
is in a state of decreasing kinetic energy and during which most centers 
of action have shifted northward. 

Several interesting temporal variations of ozone at different geo- 
graphical locations are worth mentioning. One of these locations is the 
United States. In April the 300 m*atm-cm contours extended to about 20°N 
in the longitude sector of the United States. But in May a small portion 
of this contour is noted pushing northward across 30°N over Texas. In June 
this area has enlarged to cover most of the south-central United States. 

By July a greater part of the United States has an ozone concentration of 
less than 300 m* atm -cm. This feature appears to imply that low amounts 
of total ozone are advected northward into an area centered on the 95th 
meridian. A region off the east coast of the United States, during 
April showed a north-ward indentation of 300 m. atm. cm line, and in July 
was seen to be a preferred region for a southward bulge of this contour line. 
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There was only one region in the Northern Hemisphere where the 300 
m»atm«cm contour remained north of the 30*^N latitude during all four 
months. This was the region from Central Asia to Japan. The area over 
which there was very little seasonal ozone variability at 30°N was from 
western North Africa to Iran (15°W to 60°E) . 

The 260 m. atm. cm contour in April, May, and June (280 m* atm -cm in 
July) characterizes the ozone distribution in tropical latitudes. Re- 
latively high values of 0^ appear over Africa and the eastern Atlantic 
Ocean. The departures from zonal mean values are not strong, but never- 
theless suggest the possibility that this is a preferred region (3°W to 
60°E) for interhemispherical exchange of stratospheric air masses. 

In general it appears that the western part of the Northern Hemi- 
sphere has the greatest seasonal variability of total ozone and that the 
eastern part has the least variability. 

4.213 Monthly Southern Hemisphere Variability 

There are a total of eleven active ozone surface stations in the 
entire Southern Hemisphere. As a contrast, the United States and Canada 
alone have thirteen stations. Of these eleven Southern Hemispheric sta- 
tions there are four between the equator and 30°S, five between 30°S and 
60°S, and two from 60°S to the pole. With a network as thin as this only 
a very basic understanding of the ozonosphere in the Southern Hemisphere 
exists and this has primarily come from the Australian network of five 
stations. The following maps (Figs. 4.213.1a-d) permit a first look 
at the longitudinal variations, the long-wave pattern, and the gradient 
of total ozone in the Southern Hemisphere. 

In April a 'tight' gradient is located approximately along the 100th (E) 
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Fig. 4.213.1a Same as 4.211.1a except for April 1969, Southern 
Hemisphere. 

meridian between 40. and 50°S. This gradient weakens considerably in May 
but reappears and intensifies again in June and July. In May there are 
two regions of strong meridional 0^ gradients - one over the eastern 
Atlantic and the other over the eastern Indian Ocean. The gradients 
over both these regions are somewhat less than the gradient southwest of 
Australia mentioned above. In June the strong gradient southwest of 
Australia is again evident. There is another region of tightly-spaced 
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Fig. 4.213.1b Same as 4.211.1a except for May 1969, Southern 
Hemisphere. 

contours in the south central Pacific Ocean. In July we find both of 
the features remaining in essentially the same locations. Van Loon 
(1964) reports a broad westerly wind maximum at 45-50°S from 45-100°E. 
This is in more-or-less general agreement with the ozone gradients in 
this region. The closely spaced ozone contours from 140-160°W between 
25 and 40°S would indicate a mean wind maximum is located in this region. 
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Fig. 4.213.1c Same as 4.211,1a except for June 1969, Southern 
Hemisphere. 

This is a region of the earth about which little upper atmosphere wind 
data have been published. The use of total ozone gradients to locate 
systems of this extent may be of considerable worth in the understanding 
of the general circulation characteristics of the Southern Hemisphere, 
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Fig. 4.213. Id Same as 4.211.1a except for July 1969, Southern 
Hemisphere . 

The gradient southwest of Australia is nearly as strong in April as 
it is three months later in July, the middle of the Southern Hemisphere 
winter. One might infer from this that the jet stream system over this 
region is nearly of the same velocity in the late summer as in the middle 
of winter. This would support the statements on Southern Hemispheric 
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jet streams made in Chapter 4.2. The ozone contours indicate that the 
jet stream moves slowly northward from April to July as observed from 
wind analyses. 

In the Southern Hemisphere from April to July the ozone isoline 
pattern is arranged in a pattern of four to five planetary waves. This 
is similar to the wave pattern we observed in the Northern Hemisphere. 

At least for the four months discussed in this paper, there appears to 
be four 'anchor' locations for the "ozone troughs": 100°E, 175°W, 100°W, 
and a broad area in the South Atlantic. The latter location is ill- 
defined and variable on all mean maps. 

In April we see essentially that all of the 300 m»atm. cm contour is 
south of 30°S. The only exceptions to this are small areas over south 
central Australia and in the south central Pacific. This pattern appears 
to be typical for the Southern Hemisphere in early fall. The variation 
in the oscillation of the 300 m* atm* cm contour about latitude 30°S is 
less in the Southern Hemisphere than in the Northern Hemisphere. Peculiar 
seasonal variations of total ozone are indicated over Australia. In April 
a small positive anomaly is found over this continent. In May this small 
area of maximum ozone has expanded and now covers over half of the width 
of the continent (Fig. 4.213.1b). By June the extension is complete 
and the 300 m* atm* cm contour lies north of 30°S over the entire continent. 
In July the positive anomaly shrinks slightly, but extends further west- 
ward. During this, entire period (April to July) the total ozone at 105°E 
and 30°S has shown in essence no variance. The increase of total ozone 
over Australia, as winter approaches, seems to correspond to a decrease 
of total ozone over the North American Continent as summer approaches. 


One other area with an anomalous maximum should be noted. In May 
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the 300 m* atm* cm contour has a slight 'bulge' to the north of 30 S off 
the west coast of South America (Fig. 4.213.1b). By July this maximum 
has extended to 20°S and covers all of the extreme southern part of 
South America. A final region noteworthy of mention is the western part 
of the Indian Ocean. Over this entire region an ozone minimum is re- 
flected in the 300 m* atm* cm contour. 

The 260 m*atm*cm contours in April, May, and June (280 m*atm*cm in 
July) indicate relative high total ozone concentrations in the tropical 
latitudes over Africa and the eastern Atlantic Ocean. 

Basic similarities of the distribution of total ozone in the North- 
ern and Southern Hemispheres are the following: Both hemispheres show: 

(1) a preference for a four wave ozone pattern; (2) 'tight' ozone gradients 
which, we hypothesized, are directly related to jet stream system; (3) a 
general increase (decrease) of the ozone at. a given latitude as we progress 
toward the winter (summer) season; (4) regions of relatively high ozone 
concentration in the tropics over the eastern Atlantic Ocean and the 
African Continent. 

Basic dissimilarities in the distribution of total ozone are the 
following: (1) whereas two and occasionally three regions of 'tight' 

ozone gradients were observed over the Northern Hemisphere, two and 
sometimes only one such region could be seen over the Southern Hemisphere. 
This difference may partly be based on the fact that different seasons 
in each hemisphere are being compared. This should be remedied when 
NIMBUS IV data become available for a complete year; (2) there are larger 
variations of total ozone at a given latitude in the Northern Hemisphere. 
This may also be due to the out-of-phase seasonal comparisons. 

4.3 The Longitudinal Progression of Ozone 'Waves' 
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4.31 The Variation of Total Ozone at 60°S for 17 Days in May 1969 

Figure 4.31 shows a time-longitude distribution of total ozone for 
60°S. This particular latitude was selected so that comparisons could be 
made with radiance values that were computed for the first two-thirds of 
May by Fritz (1970). 

4.311 A Comparison c f Radiance and Ozone Waves 

Fritz's co...putations were for v = 669.3 cm At this particular 
wavenumber he was measuring the radiance value of the high stratosphere. 
Ninety per cent of the flux which arrives at the satellite comes from 
above 100 mb (Fritz, 1970). This is a good frequency at which to observe 
and compare ozone variations, because it is centered in a region of the 
stratosphere between 20 and 50 mb in which a large percentage of the 
total ozone is concentrated. 

Upon examination of Figure 4.31 several features are readily appar- 
ent. Both maxima and minima of ozone progress steadily to the east. 

Fritz describes the maxima and minima of radiance values as 'waves'. 

This seems to be appropriate here as well. These ozone 'waves' move in 
the same direction as the radiance 'waves' described by Fritz (1970). 
Fritz's low (high) radiance values correspond to low (high) stratospheric 
temperatures. His locations of high and low stratospheric temperatures 
are indicated by w and c, respectively, in Figure 4.31. Note that warm 
stratospheric temperatures correspond fairly well with high concentra- 
tions of ozone. Low ozone values seem to agree with cold stratospheric 
temperatures. This correspondence has been shown from earlier balloon- 
borne ozonesonde packages (Craig, 1965; Lovill and Miller, 1968). 

4.312 Progression Rates of Eastward Moving Ozone 'Waves' 


The light lines in Figure 4.31 sloping to the right with time are 
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Fig. 4.31. Tlme-longltude graph of the total ozone at 60 S for May 
1969. Total ozone values are in m* atm. cm. The W and C symbols re- 
piresent regions of warm and cold stratospheric temperatures according 
to Fritz (1970). The straight lines represent the average tilt in 
time with longitude of the respective ozone 'waves'. 
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indications of the rate of eastward movement of the 'ozone waves’. Each 
line has a number associated with it. The average rate of longitude 
movement for each wave is indicated in Table 4.312. 

Table 4.312 


Eastward Progression Rate of Total Ozone 'Waves at 
60°S for the Period 3-19 May 1969 


Wave Number 

Progression Rate 
(Degrees of Longitude/Day) 

Average Concentration of 
Ozone Within Wave (m- atm. cm) 

1 

5,6 

380 

2 

5.7 

350 

3 

7.1 

380 

4 

13.3 

350 

5 

3.8 

360 


The variability of the longitudinal progression of these waves is 
quite large, ranging from 3.8 to 13.3 degree longitude day with an 
average movement of 7 degrees longitude day Additional studies of 
this type should help relate stratospheric ozone variations to thermal 
fluctuations, and this in turn should promote a better understanding 
of the dynamics of the stratosphere. 

4.32 Hemispherical Differences in Longitudinal Progression During June 1969 
Figures 4.32.1 to 4.32.6 represent time-longitude variations of 
total ozone at 60°S, 40°S, 20°S, the Equator, 20°N, and 40°N, respectively. 
These figures are presented in order to indicate the rapid eastward move- 
ment of djmamically active ozone maxima and minima in the winter hemisphere 
and the rather slow movement of these centers in the summer hemisphere 
and near the equator. 
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Fig. 4.32.1. Time-longitude variations of the total ozone 
(m* atm* cm) during June 1969 at 6QoS. 



Fig. 4.32.2. The same as Figure 4.32.1, except 40°S. 
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Fig. 4.32.3. 


The same as Figure 4.32.1, except 20 S. 
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Fig. 4.32.4. The same as Figure 4.32.1, except the Equator. 










-52- 


The data in Figures 4.32.1 to 4.32.6 for June 1969 represent the 
first possible analysis of longitudinal variations of total ozone on a 
fine day-by-day grid. The best data of the longitudinal distribution of 
total ozone from surface stations in the past have been at the mid-latitudes 
in the Northern Hemisphere. These data at best raised more questions 
than they answered. The data on the preceeding pages represent a solid 
data bank that extends around the earth over unpopulated a reas from the 
equator to high latitudes. 

Several facts are demonstrated quite early in these figures. Firstly, 
one can see that the speed at which the 'waves* progress eastward is 
greater in the winter hemisphere (here the Southern Hemisphere) . Secondly, 
the speed of eastward progression decreases as one approaches the lower 
latitudes in the winter hemisphere (see Figs. 4.32.1-4.32.3). Thirdly, the 
intensity of the waves, as seen from the contouring intervals, increases 
as one goes toward the pole in the winter hemisphere. Finally, in the 
equatorial region and in the Northern Hemisphere (Figs. 4.32.4-4.32.6) 
there is not a strong eastward progression of waves. Any eastward progression 
at best, is only weak and sometimes is even replaced by westward progression. 

The movements of the ozone 'waves' are of course closely related to 
winds in the upper atmosphere. Thus, one observes a strong flow from 
the west in the winter season in the Southern Hemisphere. In the summer 
in the Northern Hemisphere the flow aloft is variable and generally very 
weak. In the equatorial regions Figure 4.31.4 suggests a wind component 
from the east. The ozone 'waves' progress eastward at rates as high as 
>15° longitude/day at high latitudes in the Southern Hemisphere (see Fig. 
4.21.1). However, near the equator and at 20°N (Figs. 4.31.4-4.31.5) a 
westward movement as high as 3° longitude/day may be observed. 
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A final interesting item to mention in relation to these 'waves' 
is that as one attempts to follow them in time, one sees that they 
change intensity continuously. A good example of this phenomenon is 
seen in Figure 4.31.2 which is the fluctuation at 40°S: One may 

identify a 'wave' maximum that begins on 1 June at approximately 160°W 
and follow it to 30 June at 10°E. During this time interval the 'wave' 
progressed 170° eastward halfway around the globe. This is an average 
progression rate of 5.7° longitude/day . During the interval one can 
observe that this maximum 'wave' ranged from a high value of 380 
m* atm* cm on 13 June to a low value of 320 m* atm* cm seven days later. 

As one follows the 'wave' during these decreases and increases, it is 
also readily seen that the progression rate changes just as erratically 
as the total ozone concentration. These temporal fluctuations are 
related to the motion of dynamic stratospheric systems. The 'tighter' 
ozone gradients , which may be seen in the winter hemisphere set of 
figures (Figs. 4.32.1-4.32.3), are associated with the stronger 
stratospheric baroclinic zones (see Chapter 3) which are moving east- 
ward . 

4.4 Average Variation of Total Ozone with Latitude 

Figure 4.4 represents the monthly total ozone distribution in the 
Southern Hemisphere. London (1963), basing his analysis on the ozone 
data available through 1959, arrived at an average curve for the spring 
that extended from the equator to the North Pole. His curve is plotted 
in Figure 4.4. Notice the remarkable agreement between his total ozone 
average curve for the spring and the average IRIS curves for months 
April and May. Notice also that his summer curve is consider- 
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Fig. 4.4. The variation of total ozone with latitude, from 
80°N to 80°S for April, May, June and July 1969. 

ably lower than the curve for July from IRIS. This might indicate that 
during July 1969 there was an anomalously high ozone amount at the lati- 
tudes compared. The total number of ozone data points used to construct 
the April-July line was 1.8 x 10^. Mean total ozone is given for the 
months of April, May, June, and July as well. Immediately evident from 
Figure 4.4 are increasing ozone values from the equatorial region to the 
poles in each month. 

The largest gradient of total ozone is clearly found in the middle 
latitudes (35-50°N) in the Northern Hemisphere. The location of the 
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steepest gradient in the Southern Hemisphere is not as well defined. 

Two preferred areas are indicated there. One is in the subtropics 
from 20-25°S. The other is from approximately 30-40°S. The former 
may be influenced by the change in regression coefficients used in 
the calculations. A slope of similar steepness, perhaps also some- 
what biased by the regression coefficients, appears in the curve from 
April at 18-25°N. 

Seasonal variation of total ozone is smallest in the tropics. In 
the Northern Hemisphere this occurs at 15°N, in the Southern Hemisphere 
at 20-30°S. The latter location may be influenced to a certain extent 
by the regression coefficient changeover, but nevertheless it is felt 
that the smallest variation is in this 20-30°S latitude band. 

Another interesting feature of the figure is that the minimum 
ozone values for all latitudes for all four months investigated are 

found between 5 and 10°S. It will be very interesting to see how this 

I 

minimum migrates during the other months when Nimbus IV data are 
analyzed. 

In the Northern Hemisphere the total ozone in April reaches a 
peak at 70°N; this corresponds fairly well with London's (1963) 
winter peak at 60-65°N, but not at all with his curve for spring. 

In the Southern Hemisphere June and July have peaks at 68°S and 60°S, 
respectively. April has a secondary peak at 66°S and presumably 
another one south of 78°S. The unusual feature of an ozone minimum 
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value at 70°S is difficult to explain unless it is due to an anomalous 
standing-wave situation near Antarctica. 

Several important facts have been elaborated upon in Figure 4.4: 

(1) for the first time the average distribution of ozone with latitude 
has been determined by satellite for a period of several months; (2) 
the distribution agrees very well with the data from the long-term 
data base in the Northern Hemisphere; (3) the latitudinal ozone 
distribution, for the first time, has been described during the 
autumn and winter in the Southern Hemisphere; (4) the smallest 
variation of ozone from spring to summer was seen at 15°N; (5) the 
lowest mean ozone value for the period was at 6°S. 

4 . 5 The Global Distribution of Ozone from 19 April to 21 July 1969 

The average distribution of total ozone for the Northern and South- 
ern Hemispheres from 19 April to 21 July 1969 is given in Figures 4.5.1 
and 4.5.2. The Northern Hemisphere analysis was computer analyzed 
using 0.83 x 10^ total ozone data points. The Southern Hemisphere 
analysis used 1.01 x 10^ points. 

Several interesting features in the total ozone distribution are 
seen in these average global maps. A large region of low total ozone 
(< 240 m*atm*cm), 10 to 20 degrees on either side of the equator is 
observed to extend from southeast Asia eastward to the central Pacific. 
One region of several thousand square kilometers extent (160°E, 10°S) 
has values less than 220 m* atm* cm. Climatological surveys of this 
region in the western equatorial Pacific have shown anomalously high 




Fig. 4.5.1. The average distribution of total ozone 
(m* atm- cm) for the Northern Hemisphere from 19 April 
to 21 July 1969. Shading represents total ozone values 
less than 240 m. atm- cm. 
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amounts of low-level cloudiness and precipitation (see the work of 
Budyko (1963) on the anomalous heat balance in this area) . From earlier 
satellite studies Vonder Haar and Suomi (1971) have shown that this same 
region also has an anomalously high albedo . This would correlate well 
with the high amount of cloudiness. The manner by which this correlates 
with low total ozone concentrations is not at all clear, since most of 
the dally and seasonal change in total ozone occurs in the stratosphere 
(Craig, 1965; Lovill and Miller, 1968). There is a possibility that on 
the large-scale this region is subject to a stratospheric cell that 
allows for generally ascending motions and horizontal out-flux of 
stratospheric ozone from the area. When a complete understanding of 
the mechanism in this anomalous ozone region is reached, it will probably 
not deviate far from the above speculation, for this is obviously not a 
problem of photochemistry alone but one of atmospheric dynamics. 

A region of anomalously high total ozone in the tropics is seen to 
extend from the western Atlantic eastward across Africa to the western 
Indian Ocean. This region extends into both hemispheres. This area of 
relatively high total ozone may be associated by some interconnected 
mechanism with the anomalously low total ozone region over southeast Asia 
which is also in the tropical latitudes. 

Both the Northern and Southern Hemispheres appear to provide 
'anchoring' mechanisms for ozone ridges (high ozone amounts). The most 
predominant of these ridges is also usually associated with a 'tight' 
ozone gradient. In the Northern Hemisphere these locations are eastern 
North America, central Europe, western Asia, and eastern Asia. The 
ridge over the Soviet Union appears to be the weakest of the four ridges. 
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In the Southern Hemisphere ridges are evident over the east Indian Ocean, 
the central Pacific Ocean and the eastern Atlantic. In general the ridge 
pattern is easier to. delineate in the Southern Hemisphere than in the 
Northern Hemisphere. This may be due to a basic difference in the general 
circulation of the Southern Hemisphere or to the time period over which 
the data are analyzed. Certainly some of the difference is attributable 
to the stronger winter circulation in the Southern Hemisphere as compared 
to the weaker summer circulation in the Northern Hemisphere. 

The ’tightest’ gradient of total ozone in either hemisphere is seen 
over Japan (Fig. 4.5.1). This correlates closely with the anamalously 
high frequency of jet stream winds usually observed over the Japanese 
Islands. In the Southern Hemisphere the most tightly packed ozone contours 
are observed the eastern Indian Ocean. Another closely spaced group of 
ozone contours is seen over the Pacific Ocean, from approximately 110°W 
to 170°W, near latitude 35°S. These locations in the Southern Hemisphere 
should coincide with the highest mean winds in the upper troposphere and 
stratosphere during the period of investigation (refer to section 4.21). 

A closed center of high total ozone content (400 m*atm-cm) is 
seen in Figure 4.5.2 over the Australian sector of Antarctica. 

An average global value of the total ozone for the period from 19 
April to 21 July 1969 was calculated to be 330.3 m-atm«cm . 

4 . 6 Spectral Analyses of Ozone and Wind Variability 

Daily maps of the global distribution of total ozone were produced 
for the period from 19 April to 21 July 1969. Since only ninety-four 
days of IRIS ozone data are available at the present time, a data ensemble 
was constructued with 20° longitude intervals at various latitudes in the 
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Northern Hemisphere. In this manner a total of 1692 'unfolded days' 
of data were useable for spectral analysis. In a similar manner NMC 
250 mb wind data were analysed for the same period . 

Two examples of the spectra obtained at 50°N are presented in Figures 
4.6.1 and 4.6.2. Figure 4.6.1 shows the spectral analysis of total ozone. 

A definite spectral peak is seen at a frequency of one cycle per 12.2 
days. Figure 4.6.2 represents the spectral analysis of the 250 mb wind 
component. A spectral peak is also seen in this analysis occurring at 
one cycle per 12.2 days. A cross-spectral analysis of the total ozone 
and wind velocity was processed for several of the Northern Hemisphere 
latitudes. Table 4.6 presents these results. 


Table 4.6 

Periods at selected latitudes where the 
coherence-squared value is a maximum. 


Latitude 

Period (Peaks) 
Days 

Coherence^ 

30 N 

14.2 

0.17 

40 N 

12.2 

0.14 

50 N 

12.2 

0.36 

60 N 

43.0 

0.22 

70 N 

8.3 

0.26 ' 


An examination of Table 4.6 reveals that the 12.2 day periodicity 
at 50°N of both the total ozone and the 250 mb wind velocity has a 
coherence-squared statistic of 0.36. This indicates that the two para- 
meters are significantly related at periodicity of slightly less than 
once per two weeks. It is worthy to note that similar periodicities 
(with lesser coherence-squared values) exist also at the lower latitudes 
- 30 and 40°N (Table 4.6). This is certainly a distinct suggestion 
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Fig. 4.6.1. Spectra of the total 
ozone at 50°N from 21 April to 21 
July 1969. 


Fig. 4.6.2. The same as Figure 
4.6.1 except 250 mb wind velocity. 
The -2 line indicates the slope 
of the exponential function. 


that ozone maxima are occurring simultaneously with upper atmosphere 
wind maxima at a significant number of grid locations in the middle and 
upper latitudes of the Northern Hemisphere (at least during the late 
spring and early summer of 1969) . 





5.0 Summary and Conclusions 

Studies in recent years have shown distinct differences of the 
atmosoheric general circulation of the Northern and Southern 
Hemispheres. Earlier publications postulated a more zonal flow in 
the Southern Hemisphere supposedly because of less land mass and 
fewer mountain massifs. More recent investigations show that such 
statements need to be qualified. 

Wooldridge and Reiter (1970) have shown that significantly 
stronger horizontal anisotropy of flow prevails at cyclone wavelengths 
in the Southern Hemisphere than in the Northern Hemisphere. In the 
Southern Hemisphere meridional perturbations exceed zonal flow perturba 
tions. Eulerian spectral densities of the zonal component in the 
Southern Hemisphere were seen to be only half of those in the Northern 
Hemisphere. This would seem to imply that there is less zonal energy 
at cyclone-wave frequencies in the Southern than in the Northern 
Hemisphere. This may occur because of a lack of orographically induced 
planetary longwaves in the Southern Hemisphere that are superimposed 
upon the hemispheric jet stream zones. Apparently the relatively 
zonal character of the Southern Hemisphere pertains to the long 
planetary-scale waves only. 

Interest in global total ozone measurements is widespread and 
increasing. The paucity of total ozone measuring stations over the 
oceans and in the tropics in the Northern Hemisphere and throughout 
the Southern Hemisphere has greatly hindered a complete understanding 
of the role that stratospheric ozone plays in the general circulation 
processes of the atmosphere. With total ozone sensors (such as the 
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IRIS instrument) on earth-orbiting satellites the ozone over the 
entire earth is measured daily. The conclusions of some of the 
analysis from the first satellite (Nimbus III) capable of these 
measurements follow. 

A statistical analysis of the total ozone data for two days 
indicated that average standard error between data obtained by the 
inversion technique and data obtained by regression coefficients was 
only 3.7 percent. When surface total ozone measurements were compared 
x^ith regression data at monthly intervals a standard error of 3.6 
percent was seen. 

A strong correlation was discovered between the meridional gradient 
of total ozone and the wind velocity in jet stream systems. A graph 
has been constructed which will allow the computation of the wind 
velocity at 250 mb if the total ozone distribution in the vicinity is 
known. This should be a very useful relation to apply over the oceanic 
regions of the Northern Hemisphere and throughout the entire Southern 
Hemisphere. Wind velocities at jet aircraft operation levels could 
then be forecasted by consulting the daily satellite total ozone 
analysis . 

A study of the total ozone distribution over a Pacific Ocean extra- 
tropical cyclone and its associated jet stream indicated large vertical 
motions and strong horizontal advective processes were present near this 
system. Large amounts of total ozone were measured immediately behind 
the system and low amounts in front which were indicative of air rich in 
ozone moving southward and downward behind the system and air poor in 
ozone moving northward and upward in front of the storm. 
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Time-longitude variations of the total ozone for several 1 ititudes 
in the Northern and Southern Hemisphere for June 1969 were analyzed. 
Several interesting facts were seen in these results. Firstly, one can 
observe that the speed at which the 'waves' progress eastward is greater 
in the winter hemisphere (here the Southern Hemisphere). Secondly, the 
speed of eastward progression decreases as one approaches the lower 
latitudes in the winter hemisphere. Thirdly, in the equatorial region 
and in the Northern Hemisphere there is not a strong eastward progres- 
sion of 'waves' but at best only weak eastward movement. Progression 
rates as high as fifteen degrees of longitude per day were seen at the 
high latitudes in the Southern Hemisphere. As one follows the 'waves' 
in time, continuous changes in intensity are observed. 

When total ozone data for the period of observation are plotted 
against latitude for the various months the following is indicated. 

For the first time the average distribution of ozone has been described 
during the autumn and winter at close intervals, both temporally and 
spatially, in the Southern Hemisphere. The smallest variation of 
total ozone from spring to summer was seen at 15°N. The lowest mean 
total ozone value for the globe was at 6°S. 

When two maps were prepared describing the distribution of total 
ozone for the entire period for the Northern and Southern Hemispheres, 
the following was seen. Closely spaced total ozone isolines appear to 
have revealed a reasonably strong jet stream west of South America at 
approximately 35°S. This system, at least at the intensity indicated 
by the total ozone gradient, has gone more-or-less unnoticed in the 


past . 



-66 


Additionally, a large region of low total ozone was observed In 
the tropics, extending from southeast Asia eastward to the central 
Pacific. Earlier studies have Indicated that this region also has an 
anomalously high albedo. The manner by which this correlates with low 
total ozone concentrations is not at all clear . 

Both the Northern and Southern Hemispheres appear to provide 'an- 
choring' mechanisms for total ozone ridges. The most predominant of 
the ridges is also usually associated with a 'tight' ozone gradient. 

In the Northern Hemisphere these locations are eastern North America, 
central Europe, western Asia, and eastern Asia. The ridge over the 
Soviet Union appears to be the weakest of the four ridges. In the 
Southern Hemisphere 'ozone' ridges are evident over the east Indian 
Ocean, the central Pacific Ocean and the eastern Atlantic. In general, 
the ridge pattern is easier to delineate in the Southern Hemisphere than 
in the Northern Hemisphere . 

The 'tightest' gradient of total ozone in either hemisphere was 
seen over Japan. This correlates closely with the high frequency of 
jet streams observed over the region. 

An average global value of the total ozone for the entire period 
3JO.J 

was calculated to be MKtKK m*atm<cm. 



-67- 


5.1 Suggestions for Future Research 

The results presented In this paper represent a significant advance 
of our knowledge toward a better understanding of the manner in which 
ozone is involved in depicting both the large- and small-scale motions 
of the atmosphere. Future work in this area should take advantage of 
the ever-improving instrumentation on current and future satellites in 
order to improve global coverage. Many of the results can be significantly 
amplified in detail when the Nimbus IV IRIS data are released. All four 
seasons in both hemispheres will then be available for analysis and 
comparison. We should be able to determine accurately variations in 
location and strength of major motion systems in the atmosphere, such 
as jet streams, cyclones, and anticyclones. More proof will be offered, 
for example, to investigators who have suggested that stronger horizontal 
anisotropy of flow prevails at cyclone wavelengths in the Southern 
Hemisphere rather than in the Northern Hemisphere. The techniques 
developed in this paper should be useful for the analysis of data not 
only from earth-orbiting satellites but from those orbiting other planets, 
such as the current Martian probe by the Mariner 9 IRIS experiment. 
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APPENDIX II 


Continued 


Documentation o£ Major Computer Programs 
Developed under grant auspices 

Program TIMELON 

Program 03C0NT 



PROGRAM: TIMELON 


This program contours an array formed by taking the ozone values 
along a specified latitude line on succeeding days and placing the 
"lines" in the array by their day number . The analysed data is 
plotted on microfilm with longitude on the horizontal axis and time 
on the vertical axis. The time period covered can be varied in units 
of 1 day up to the total period for which data is available. Each 
frame of microfilm contains 1 labeled graph. 

Input Cards: 

Card 1 


Description 


Columns 


1-4 

LDAY 

5-8 

KDAY 

Card 2 and 3 


Columns 


1-80 

with format 2014 

NODAY 

Card 4 


Columns 


1-80 


Card 5 



Integer, the day no. of the first day 
to be included in the graph. 

Integer, the day no. of the last day 
to be included in the graph. 


Integer, an array of 23 elements with 
the day nos. of the days that have no 
data. 


Descri pti on 

Alphanumeric label for microfilm 


Columns 

1-50 


Columns 

1-20 


Description 
more labels 


Description 
more labels 


Description 


Card 6 



Method 


The ozone values are read from the data tape along with their 
associated day no.* latitude and longitude. Observations outside 
the time period of the graph to be generated or with a latitude 
greater than 15°from the desired latitude line and had data are 
eliminated. The ozone values for the remaining observations are 
fitted to a Rat-Ion grid with 2° increments. After all the obser- 
vations for one day are processed the grid is interpolated and the 
values along the specified latitude line are placed in another array 
according to the day number. The frid is then zeroed and the next 
days observations are processed. When the observations for the 
entire time period have been processed, the second array is inter- 
polated to supply the values for (possible) missing days. The array 
is contoured and the results are plotted on microfilm. 



PROGRAM 03C0NT 


This program objectively contours ozone data and outputs the 
analysed data on microfilm. The data is contoured over a world 
map (cylindrical equidistant projection) showing continental outliner 
and with latitude and longitude lines labelled. Each frame of microfilm 
contains one world map. The time period covered by each map can be 
verified in units of 1 day from 1 day to the total period for which 
data is available. 


Input Cards: 

Card 1 

Columns 

(L)line 

1-4 


5-8 


9-14 


two blank lines 
Card 2 and 3 
blank line 
Col urns 
1-80 


Card 4 
Columns 
1-80 
Card 5 
Columns 
1-50 
Card 6 
Columns 


Description 

LDAY Integer, day no. of the first day of the 
time period to be covered by the map. 

KDAY Integer, day no. of the last day of the 
time period to be covered by the map 
(may be same as LDAY). 

LENGTH Integer, length of the array to be read 
from take containing plotter instructions 
for generating the background map. 


Descri pti on 

NODAY Integer, an array of 23 elements containing 
the day nos. of the days that have no data. 


Description 

A Yehanumeric label for microfilm. 


Descri pti on 
more labels 


1-20 


Descri pti on 
more labels 



Method 


The ozone values are read from the data take along with their 
associated day no latitude and longitude. The day no. is checked 
to determine if the observation is during the period under consider- 
ation. Then the ozone value is tested against extreme reasonable 
values and had data is eliminated. The observations passing the 
above tests are fitted to a Rat-Ion grid with 20 increments between 
grid lines. When all the observations in have been processed, the grid 
is interpolated linearly by rows and by columns. The results are then 
merged weighting the rows 10 to 1 over columns. The grid is then con- 
toured and the results plotted on microfilm. 



